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ABSTRACT 


A  reinvestigation  was  conducted  of  an  anomaly  in  the  experimental 
angular  correlation  for  the  gamma-gamma  cascade  of  Ni^O  which  was 
reported  by  Verser  and  Conway  in  1960  and  1961,  respectively.  The  measure¬ 
ment  of  this  correlation  was  being  used  to  test  angular  correlation 
equipment  built  at  the  US  Naval  Postgraduate  School.  Two  new  experimental 
determinations  of  the  correlation  function  with  the  same  equipment  as  used 
previously  indicated  good  agreement  with  the  theoretical  function. 

Verser' s  and  Conway's  experimentally-determined  angular  correlation 
functions  were  found  to  deviate  less  significantly  from  the  theoretical 
function  than  they  had  reported  since  their  standard  deviations  had  failed 
to  take  into  account  the  variations  of  the  characteristics  of  the  electronic 
equipment  due  solely  to  electronic  causes.  It  is  concluded  that  the  anomaly 
reported  by  previous  experimenters  was  statistical  in  nature. 

Programs  for  the  CDC  1604  computer  have  been  developed  in  FORTRAN 
language  for  complete  data  reduction  accepting  as  input  the  total  counts 
in  each  channel,  the  total  coincident  counts,  and  other  experimental 
parameters.  The  program  provides  for  the  calculation  of  the  coefficients 
for  a  cosine  series  and  a  Legendre  polynomial  series  expansion  of  the 
least  squares  fit  to  the  experimental  points.  Other  programs  were 
developed  to  calculate  the  corrections  for  non-point  detectors  to  be 
applied  to  the  theoretical  coefficients,  and  to  calculate  coincidence 
resolving  times  from  data  obtained  during  a  large  number  of  runs. 

Major  improvements  were  made  in  the  equipment  to  be  used  in  future 
angular  correlation  experiments,  including  the  construction  and  testing 
of  both  a  slow  and  a  fast  coincidence  unit. 

The  author  wishes  to  express  his  appreciation  to  Associate  Professor 
H.  E.  Handler  who  has  guided  the  development  of  angular  correlation 
measuring  equipment  at  the  US  Naval  Postgraduate  School  for  the  last 
three  academic  years.  His  recommendations,  particularly  in  the  area 
of  data  analysis,  and  his  perusal  of  the  draft  of  this  paper  resulted  in 
a  still  more  thorough  reinvestigation  of  the  results  of  the  angular 
correlation  studies  reported  by  previous  experimenters  at  this  School. 

It  is  a  pleasure  to  ackowledge  the  guidance  and  aid  received  from  Lt  Cmdr 
P.  A.  Phelps  (USN)  in  comprehending  the  theoretical  aspects  of  the  basic 
experiment,  in  becoming  familiar  with  the  electronic  equipment,  and 
in  writing  and  debugging  the  computer  programs.  Appreciation  is  expressed 
to  Mr.  M.  A.  Brillhart  for  his  construction  of  the  coincidence  units,  and 
for  his  assistance  in  developing  the  mixer  circuitry. 
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INTRODUCTION 

A.  Purpose  of  Angular  Correlation  Experiments 

In  gamma-gamma  angular  correlation  experiments,  one  observes 
the  probability  of  emission  of  the  second  quantum  of  a  cascade  as 
a  function  of  angle  relative  to  the  first  quantum.  In  general, 
this  probability  is  an  anisotropic  function  of  the  angle.  Such 
experiments,  in  conjunction  with  other  nuclear  studies,  can  lead 
to  identification  of  the  multipole  orders  of  the  emissions  and 
the  total  angular  momenta  of  the  original,  intermediate,  and  final 
states  of  the  nuclide  under  investigation.  ^1,2,3  The  orienta¬ 
tion  of  the  total  angular  momentum  vector  of  the  intermediate  level 
is  correlated  with  the  direction  of  emission  of  each  quantum.  The 
anisotropy  of  the  probability  distribution  pattern  will  be  attenuated 
if  the  orientation  of  the  total  angular  momentum  vector  of  the  inter¬ 
mediate  level  changes  between  the  times  of  these  two  emissions.  Thus, 
one  must  study  a  nuclide  for  which  the  intermediate  energy  level  has 
a  mean  life  which  is  short  relative  to  the  period  of  precession  of 
the  angular  momentum  vector  of  that  level.  Only  if  the  total  angu¬ 
lar  momentum  vector  of  the  intermediate  level  does  not  appreciably 
change  its  orientation  between  the  emissions  of  the  two  quanta  can 
the  probability  distribution  as  a  function  of  angle  between  the  two 
emissions  be  interpreted  meaningfully. 

Generally,  no  parity  assignments  can  be  made  for  the  nuclear 
levels  involved  based  solely  on  a  directional  angular  correlation 
study  such  as  described  above,  since  the  angular  distribution  of 
the  emissions  is  the  same  for  electric  and  magnetic  multipole  radi¬ 
ations  of  the  same  order.  By  means  of  a  similar  polarization- 
direction  correlation  experiment,  wherein  one  measures  the  orientation 
of  the  polarization  vector  relative  to  the  plane  of  the  two 

quanta,  one  can  determine  whether  the  emissions  are  electric  or  mag¬ 
netic  multipoles.  [ 4 , 5 j  This  information,  coupled  with  known  mul¬ 

tipole  orders  obtained  from  an  angular  correlation  study,  enables 
the  assignment  of  relative  parities  to  the  three  nuclear  energy 
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levels.  Absolute  parity  assignments  require  additional  data* 

In  the  case  of  Ni^,  an  even-even  nuclide,  an  even  parity 
and  zero  total  angular  momentum  are  predicted  for  the  ground 
state  on  the  basis  of  the  nuclear  shell  model,  and  are  corrobo¬ 
rated  by  other  experimental  techniques  unrelated  to  this  paper . 
Studies  of  the  decay  scheme  have  shown  that  the  spins  of  the 
first  and  second  excited  levels  are  2  and  4  respectively .  Cor¬ 
relation  studies  are  consistent  with  this  and  demonstrate  that 

these  radiations  are  pure  electric  quadrupoles.  Thus  the  com- 

60 

plete  decay  scheme  for  Ni  is 


Decay  Scheme  of  Ni 
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Since  this  decay  scheme  is  well  understood,  this  gamma-gamma 
cascade  of  Ni^  is  frequently  used  to  determine  whether  angular 
correlation  equipment  is  functioning  properly*,  f  6  j  It  is  used 
in  this  manner  in  this  paper. 

B.  Methods  of  Conducting  the  Experiment 

Since  nuclei  are  generally  randomly  oriented,  one  cannot 
measure  the  radiation  pattern  of  a  given  transition.  One  way  of 
minimizing  the  random  orientation  of  the  nuclei  is  to  apply  a 
strong  magnetic  field  or  a  strong  electric  field  gradient  at  very 
low  temperatures;  radiation  patterns  have  been  measured  using  this 
technique.  In  the  event  of  a  cascade,  the  effects  of  this  random 
orientation  can  be  partially  circumvented  by  establishing  some  re¬ 
ference  direction  with  respect  to  which  only  certain  of  the  second 
of  the  two  radiations  are  measured.  In  the  experiment  discussed  in 
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this  paper,  nuclei  have  not  been  aligned;  the  reference  direction 

for  observation  is  the  direction  of  emission  of  the  first  quantum 

60 

in  the  two-gamma  cascade  of  Ni  ,  Since  the  lifetime  of  the  intern- 
mediate  nuclear  level  for  this  nuclide  is  short  in  comparison  with 
the  period  of  precession  of  its  angular  momentum  vector,  the  ori- 
entation  of  the  angular  momentum  vector  does  not  vary  significantly 
between  the  times  of  emission  of  the  two  quanta •  Consequently,  an 
unattenuated  angular  correlation  between  the  directions  of  emission 
of  the  two  successive  radiations  is  observed. 

C.  Interpretation  of  the  Experimental  Data 

By  measuring  the  total , coincidence  counting  rate  at  several 
angles  9. between  successive  quanta,  and  applying  corrections  for 
accidental  coincidences  and  for  the  attenuation  of  the  correlation 
function  due  to  the  use  of  detectors  of  finite  size,  7,8*  one  can 
determine  the  relative  likelihood  of  emission  of  the  second  quantum 
at  these  angles  relative  to  the  direction  of  emission  of  the  first  a 
The  experimental  estimate  of  the  correlation  function  at  a  specific 
angle  0^  is  calculated  according  to  the  formula 


W  (6.) 
ex  i' 


N1N2 


(1) 


where  C  is  the  genuine  coincidence  counting  rate  at  the  given  angle, 
G 

and  Nj  and  ^  are  the  individual  channel  counting  rates  at  that 

angle  corrected  for  the  system  dead  time  and  background.  The  least 

squares  fit  to  these  data  is  the  experimental  correlation  function 

Wis(9).  7  W^s(0)  exPressed  in  terms  of  a  linear  combination 

of  even  powers  of  cos  9  or  in  terms  of  a  linear  combination  of  even 

order  Legendre  polynomials.  In  the  case  of  the  gamma  gamma  cascade 
60 

of  Ni  ,  the  theoretical  correlation  function,  W^^CO) 3  corrected  for 

the  solid  angles  subtended  by  the  detectors  used  in  this  experiment 

* 

is  given  by 


wth(e>  =  1.00000  +  0.09960  P2(cos  0)+  0.00838  P4(cos  0)  .  (2) 
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Experimental  results  can  be  compared  with  theory  by  comparing  the 
coefficients  calculated  in  the  least  squares  fit  with  the  above 
coefficients.  In  addition,  one  can  compare  the  experimentally- 
determined  anisotropy  with  the  theoretical  value.  The  anisotropy 
is  defined  as  follows. 


Wls(180°) 
“  Wls(90°) 


(3) 


The  theoretical  value  of  R,  corrected  for  the  geometry  of  this 
experiment  is  +  0. 16221 .  A  close  agreement  between  the  experi¬ 
mental  and  theoretical  coefficients  and  between  the  experimental 
and  theoretical  anisotropies  provides  a  check  of  proper  equipment 
operation* 

Do  Previous  Work  in  This  Field  at  U.So  Naval  Postgraduate  School 
The  basic  equipment  for  this  experiment  was  assembled  in 
academic  year  1959-60  by  Verser  [  9  ,  who  conducted  an  angular 
correlation  experiment  in  the  90-180°  quadrant  and  reported  agree¬ 
ment  with  the  theoretical  anisotropy  but  reported  marked  deviation 
from  theory  in  the  shape  of  the  experimental  curve.  The  equipment 
was  modified  slightly  in  academic  year  1960-61  by  Conway  £loj  who 
performed  the  same  experiment  in  the  same  quadrant  and  also  reported 
a  similar  deviation  from  theory  in  the  shape  of  the  experimental 
curve.  Fig  (2)  shows  the  results  of  these  investigations.  Verser 
reported  a  dip  in  his  experimental  curve  which  deviated  by  more  than 
two  standard  deviations  from  the  theoretical  curve  in  the  120-135° 
region.  Conway  reported  a  similar  dip  in  the  vicinity  of  135°,  and 
a  less  pronounced  dip  in  the  120°  region, 

Conway  wrote  two  programs  in  CDC  1604  computer  assembly  routine 
language,  one  for  the  least  squares  fit  of  a  linear  combination  of 
Legendre  polynomials  of  even  order  to  the  experimental  data,  and  a 
second  to  calculate  the  solid  angle  corrections  necessitated  by  the 
finite  size  of  the  detectors.  The  least  squares  program  required  a 
considerable  amount  of  hand-calculated  data  reduction  before  it  could 
be  used. 
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Possible  mechanisms  underlying  the  reported  anomalous  dip 
in  the  experimental  curve  as  suggested  by  Verser  and  Conway  in~ 
elude  the  following: 

1,  A  change  in  the  conditicnof  the  source,  possibly 
crystallization  resulting  in  preferential  align¬ 
ment  of  the  nuclei  within  the  crystal  due  to  in¬ 
ternal  electric  fields . 

2,  Fluctuations  of  ambient  temperature  which  resulted 
in  changes  in  the  operating  characteristics  of  the 
equipment „ 

3,  Uncompensated  variations  in  line  voltage, 

4,  Perturbing  effects  of  external  magnetic  fields, 

5,  Improper  orientation  of  the  source  relative  to  the 
detectors . 

It  was  recognized  that  the  likelihood  of  having  any  noticeable 
effects  due  to  Items  "l"  or  ,,4"  above  at  room  temperature  was  ex¬ 
tremely  small  since  the  Boltzmann  factor  for  the  relative  popula¬ 
tions  of  the  energy  levels  is  nearly  unity  except  under  the  influ¬ 
ence  of  extremely  strong  fields. 

In  addition,  the  possibility  was  recognized  of  there  being  in 
fact  no  real  anomaly,  but  that  the  discrepancy  was  simply  an  un¬ 
likely  statistical  occurrence  or  that  some  systematic  error  existed 
in  the  experimental  technique,  in  the  processing  of  data  by  hand,  or 
in  one  of  the  previous  computer  programs.  The  likelihood  of  this 
deviation  from  the  theoretical  being  simply  a  statistical  fluke  be¬ 
came  greater  when  it  was  recognized  that  the  variation  in  the  re¬ 
solving  time  of  the  coincidence  unit  was  not  taken  into  account  in 
computation  of  the  standard  deviation  of  the  correlation  function 
evaluated  at  the  angles  at  which  the  experiments  were  conducted, 

E,  Scope  of  the  Investigation 

Many  experimenters  have  reported  agreement  with  the  theoretical 
gamma-gamma  angular  correlation  function  for  Ni^  j^ll,12j,  yet  two 
experimenters  at  this  School  independently  reported  a  significant 
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deviation  from  theory  in  the  105-135°  region.  Thus  it  was  clear 
that  before  one  could  use  this  equipment  for  the  experiments  for 
which  it  was  assembled;  namely,  experiments  with  short-lived  nu¬ 
clides,  further  investigation  was  necessary  either  to  arrive  at 
a  mechanism  for  the  reported  anomaly,  or  to  demonstrate  that  the 
anomaly,  under  more  closely  controlled  experimental  conditions, 

is  not  observed.  In  the  investigation  of  this  anomaly,  certain 

1 

improvements  were  made  in  the  experimental  equipment;  new  pro- 

i 

grams  in  FORTRAN  language  were  developed  to  compute  the  solid 
angle  corrections  to  the  angular  correlation  function  and  cal¬ 
culate  the  least  squares  fit;  and  experiments  were  conducted 
under  closely  controlled  conditions  in  both  the  90-180°  and  the 
180-270°  quadrants  to  see  whether  a  similar  dip  would  be  observed. 

Two  additional  coincidence  units  were  constructed.  One  of 
these  is  a  slow  coincidence  unit  which,  unlike  that  used  by  Verser 
and  Conway,  has  a  resolving  time  which  is  independent  of  counting 
rate.  This  characteristic  will  be  required  in  conjunction  with 
future  experiments  with  short-lived  nuclides.  The  second  coinci¬ 
dence  unit  has  a  resolving  time  an  order  of  magnitude  shorter  than 
that  of  the  slow  units.  The  use  of  the  fast  unit  in  conjunction 
with  a  slow  unit  enables  one  to  achieve  better  precision  by  reducing 
the  accidental  coincidence  counting  rate. 
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II 


REINVESTIGATION  OF  THE  EARLIER  REPORTED  ANOMALY 

A.  Equipment  Used 

The  experimental  setup  employed  by  Verser  and  Conway  is  des¬ 
cribed  in  detail  elsewhere  |^9,loj.  Figure  (3)  indicates  the  general 
test  setup  and  Figure  (4)  shows  the  angular  correlation  table  with 
the  source,  crystals,  photomultiplier  tubes,  and  preamplifiers.  The 
same  slow  coincidence  unit  as  used  by  Verser  and  Conway  was  used  in 
the  angular  correlation  experiments  in  the  90-180°  and  in  the  180- 
270°  quadrants.  The  only  significant  change  in  the  setup  for  these 
experiments  consisted  of  adding  a  voltage  regulator  to  minimize  the 
effects  of  uncompensated  line  voltage  fluctuations .  Regulated  vol¬ 
tage  was  then  supplied  to  those  parts  of  the  electronic  apparatus 
most  likely  to  distort  the  experimental  measurements  when  under  the 
influence  of  line  voltage  fluctuations.  These  included  the  coinci¬ 
dence  unit,  both  linear  amplifiers,  and  the  high  voltage  power  supply 
for  the  photomultiplier  tubes. 

B.  Coincidence  Resolving  Time  Tests 

The  accidental  coincidence  counting  rate  is  calculated  from  the 
usual  formula 


CA  =  2  r  NjN2,  (4) 

where ^  is  the  coincidence  resolving  time.  Before  and  during  the 
angular  correlation  experiments,  scores  of  tests  were  conducted  to 
determine  the  stability  and  the  magnitude  of  T  .  The  unit  exhibited 
random  fluctuations  in  resolving  time  as  much  as  three  times  as  large 
as  those  that  can  be  explained  in  terms  of  nuclear  statistics.  The 

r  ClD 

fractional  variation  in  '  ,  y  was  observed  to  be  approximately 
twice  that  reported  by  Verser  and  Conway.  The  measured  mean  is 
0.2470  microseconds.  The  variance  of  a  single  measurement  of  T  is 
0.0044  microseconds;  this  includes  the  variation  due  to  non-nuclear 
fluctuations  and  was  used  in  calculating  the  variance  of  the  acci¬ 
dental  coincidence  counting  rates  in  this  work.  Since  the  resolving- 
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Figure  3,  Experimental  Setup  of  Verser  and  Conway 


Xtal*  s 
PM' s 

Amp1  s 


Coinc 


=  Harshaw  Scintillator  Crystals,  Type  8S8/2,(2n  diam  x2n  thick) 

=  RCA  6342  Photomultipler  Tubes 

Power  for  PM  tubes  from  Hamner  H.V\  Power  Supply  Model  N-401 

-  Hamner  N302  Amplifiers  -  each  includes  a  differential  pulse 
height  analyzer  which  supplies  a  standard  pulse  out  ("PHS"  out) 
(A  linearly  amplified  input  is  also  supplied  as  an  output  but 
was  not  used  in  this  experimental  setup) . 

=  Coincidence  Analyzer  *  Atomic  Model  502A. 


Scalers  -  Baird-Atomic,  Model  1283 
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Fig.  4  Angular  Correlation  Table 


time  characteristics  were  investigated  and  documented  far  more  com 
pletely  than  in  previous  years,  the  observed  instability  was  felt 
to  be  similar  to  that  which  would  have  been  found  by  the  earlier 
experimenters  had  they  investigated  this  problem  more  thoroughly. 

An  inclusion  of  the  effect  of  this  fluctuation  in  the  analysis  of 
their  data  would  decrease  their  reported  precision. 

C.  Daily  Checkout  Procedure  and  Instrument  Log 

A  daily  checkout  procedure  was  devised  to  ensure  that  exactly 
the  same  checks  of  the  operating  characteristics  of  the  electronic 
equipment  were  performed  each  day  in  the  same  sequence  and  with  the 
equipment  in  the  same  operating  configuration  at  each  stage  of  the 
checkout . 

v  An  instrument  log  was  developed  in  which  the  experimental 
equipment  layout  was  recorded,  and  monitoring  points  for  key  wave 
shapes  were  indicated.  Photographs  of  these  key  wave  shapes  are 
included,  as  well  as  the  record  of  the  oscilloscope  settings  used 
when  the  photographs  were  taken;  a  sample  page  from  this  log  is 
shown  in  Figure  5  °  This  log  has  proven  useful  in  trouble-shooting 
the  extensive  circuitry  and  for  monitoring  changes  in  operating 
characteristics,  such  as  pulse  rise  times,  which  could  degrade  system 
performance . 

D.  Scaler  Dead  Time 

As  a  part  of  the  complete  checkout  and  alignment  of  all  elec¬ 
tronic  gear  prior  to  conducting  the  angular  correlation  experiments, 
two-source  experiments  were  run  to  determine  the  system  resolving 
time  (primarily  scaler  dead  time) .  Several  of  these  measurements  all 
led  to  the  same  surprising  result,  namely,  a  negative  dead  time  for 
the  system.  Further  investigation  of  this  problem  involved  supplying 
pairs  of  pulses  from  a  double  pulse  generator  at  a  constant  repetition 
rate  and  varying  the  time  between  the  pulses  of  each  pair.  Expected 
behavior  was  that  with  the  1000  per  second  repetition  rate  for  the  pair 
of  pulses,  2000  counts  per  second  would  be  registered  on  the  scalers  as 
long  as  the  separation  between  the  pulses  of  each  pair  exceeded  the 
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resolving  time  of  the  system.  When  the  separation  between  the  two 
pulses  was  less  than  the  resolving  time  of  the  system,  1000  counts 
per  second  would  be  expected.  What  was  observed  was  the  recording 
of  2000  counts  per  second  when  the  interval  between  the  pulses  of 
each  pair  was  greater  than  about  5-10  microseconds,  then  as  the  in¬ 
terval  was  decreased,  the  count  rate  recorded  by  the  scalers  in¬ 
creased,  on  occasion  by  as  much  as  257e;  i.e.  ,  to  about  2500  counts 
per  second. 

Since  the  experimentally-determined  correlation  function  is 
not  sensitive  to  the  introduction  of  these  spurious  counts  because 
the  error  introduced  into  the  single  channel  counting  rates  is  es¬ 
sentially  constant  for  all  angles  and  the  effect  on  the  coincidence 
counting  rate  is  negligible,  this  anomalous  behavior  had  no  signifi¬ 
cant  influence  on  the  experimental  results  reported  herein.  New 
scalers  will  be  employed  for  future  experiments. 


E. 


Calculation  of  the  Standard  Deviation  of  W  (0.) 

ex  i/ 

The  method  of  calculating  the  standard  deviation  of  the  experi¬ 
mental  estimate  of  the  correlation  function  has  been  the  subject  of 
extensive  consideration.  The  experimental  estimate  of  the  correlation 
function  at  a  specific  angle  9^  is  calculated  according  to  the  formula 


W  (0.)  = 
exN 


N1N2 


(1) 


where  C  is  the  genuine  coincidence  counting  rate  at  the  given  angle, 
G 

and  and  are  the  individual  channel  counting  rates  at  that 
angle  corrected  for  the  system  dead  time  and  background.  The  genuine 
coincidence  counting  rate  is  calculated  from 


where  C  =  the  genuine  coincidence  counting  rate,  C_  -  the  total  coinci- 
G  o  T 

dence  counting  rate,  and  =  the  accidental  coincidence  counting  rate. 
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CT  and  the  single  channel  counts  were  measured  directly  by  scalers 
and  the  accidental  coincidence  rate  was  calculated  from 


CA  =  2  T  N1  N2  •  (6) 

It  can  be  seen  that  the  standard  deviation  of  W  (9.)  depends 

ex'  l 

almost  entirely  upon  that  of  C  .  The  standard  deviation  of  C  is 

G  G 

obtained  by  taking  the  square  root  of  the  sum  of  the  variances  of 
the  two  presumably  independent  quantities,  and 


(T 


(CG>  “  \ 


<r2(cT)  +<t2(ca) 


(7) 


Verser  and  Conway  calculated  this  standard  deviation  assum¬ 
ing  that  the  distributions  of  C  and  C  are  Poisson  in  nature  and 

A  X 

therefore  that  the  variances  are  expressible  as: 

(T2(Ca)=^  and  (T*<Ct)  =  (8a,  b) 

where  t  is  the  duration  of  the  measurement.  Therefore  they  calcu¬ 
lated  (C  )  from  the  relation 
G 


r(6G)  - 


(9) 


If  the  equipment  used  in  the  experiment  were  perfectly  stable,  with 

no  changes  in  T  due  to  changes  in  electronic  characteristics,  then 

all  fluctuations  in  C  and  C  would  be  consistent  with  Poisson  dis- 

A  G 

tributions  and  these  equations  would  apply.  As  pointed  out  in 
2  ’ 

Section  B,  ff~  (C )  when  calculated  from  equation  8a  is  too  small. 

To  calculate  a  realistic  variance  of  C^,  one  must  provide  not  only 
for  the  statistical  variations  due  to  the  random  nature  of  the  radio¬ 
active  decay  process,  but  also  for  the  variations  in  the  resolving 
time  of  the  coincidence  unit  due  to  changes  in  its  electronic  oper¬ 
ating  characteristics. 
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If  electronic  instabilities  were  present  and  nuclear  fluc¬ 
tuations  were  small  in  comparison  thereto,  then  the  variance  of 
would  be  described  adequately  by  the  relation 


cr) 


-i£ 


r 


b 


3 


uat  I  O  W  (  (o) 


Then 


( /  o) 


where  (T~  _  denotes  the  standard  deviation  due  solely  to  the  elec- 
elec 

tronic  instability  of  the  coincidence  unit. 

Since  under  experimental  condit ionsneither  of  the  two  limiting 

conditions  obtains;  i.e.,  perfect  equipment  stability  in  which  case 

the  variance  of  would  be  Poisson,  or  infinite  counting  time  in 

which  case  nuclear  fluctuations  would  be  insignificant  in  comparison 

with  (P  (C  ),  both  statistical  deviations  must  be  taken  into 

elec  A/  2  • 

account  in  calculating  (C  ) .  While  the  technique  for  the  rigorous 

A 

statistical  combination  of  these  variances  is  not  certain  in  spite 
of  efforts  to  arrive  at  a  valid  analytical  expression,  the  rather 
naive  approach  of  adding  them  would  seem  to  be  a  conservative  ap¬ 
proach  which  applies  to  the  experimental  situation.  Thus  the  ex- 

2  • 

pression  used  in  this  report  for  calculating  (p  (C  )  is: 

G 


It  can  be  seen  that  this  expression  gives  the  correct  standard 
deviation  in  the  two  limiting  cases  discussed  earlier,  and  it  is 
assumed  that  the  error  resulting  from  applying  this  expression  in 
the  actual  situation  is  small. 
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On  the  basis  of  the  experimental  determination  of  T  ,  the  third 
term  of  this  expression  is  of  the  same  order  as  the  first  two;  hence, 
failing  to  take  it  into  account  results  in  unrealistically  small  esti¬ 
mates  of  the  uncertainties  in  the  experiment. 

F  0  Preparation  and  Check  of  New  Computer  Programs 

Conway  prepared  two  programs  in  1960-61  in  1604  computer  assem¬ 
bly  routine  language.  The  purpose  of  one  was  to  calculate  the  solid 
angle  correction  following  the  scheme  of  Rose  j^8 j  to  apply  to  the 
theoretical  correlation  function  to  account  for  the  attenuation  re- 
suiting  from  the  finite  dimensions  of  the  detectors.  The  second 
program  calculated  the  coefficients  for  a  least  squares  fit  £ 7 j  of 

even  order  Legendre  polynomials  to  W  (0.). 

ex  l 

As  mentioned  previously,  difficulties  were  experienced  in  at¬ 
tempting  to  run  these  programs.  In  addition,  the  least  squares  pro¬ 
gram  required  the  following  prior  manual  computations 

1.  application  of  scaler  dead  time  corrections 

2.  correction  of  counting  rates  for  background 

2  • 

3.  calculation  of  CA  and  (T  (C4) 

A  A 

.  2  9 

4.  calculation  of  C  and  (p  (C  ) 

G  G 

r\  _ 

5.  calculation  of  W  (9.)  and  <p  [w  (Q.)| 

ex'  V  i  ex  1  J 

6„  calculation  of  a  weighted  mean  of  ^ex(®^)  at  each  angle 

7.  calculation  of  the  standard  deviation  of  the  weighted 
mean  at  each  angle 

The  values  of  the  weighted  mean  of  W^CO^)  and  the  standard  deviation 
of  the  weighted  mean  at  each  angle  were  the  entering  arguments  for 
Conway's  least  squares  program. 

The  original  intent  was  to  program  the  above  listed  preliminary 
computations  only,  and  to  append  Conway's  least  squares  program  as  a 
subroutine.  It  was  not  possible  to  append  his  program  to  the  FORTRAN 
program,  however,  due  to  the  incompatibility  of  the  two  computer 
languages.  A  FORTRAN  program  was  written  to  accept  as  input  the  total 
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counts  registered  by  the  scalers  and  thus  save  considerable  time 
and  achieve  greater  accuracy  by  avoiding  the  errors  likely  to  creep 
into  the  lengthy  manual  computations. 

Separate  programs  were  written  which  would  compute  the  coeffici¬ 
ents  for  the  even  power  cosine  series  as  well  as  the  coefficients  for 
the  even  order  Legendre  polynomial  series  representation  of  W^(9), 
calculate  the  solid  angle  corrections  to  be  applied  to  W^(9) ,  and 
calculate T  from  the  data  obtained  during  coincidence  resolving 
time  runs. 

The  above  programs  provided  independent  checks  of  Conway's 
least  squares  and  solid  angle  correction  programs.  The  new  least 
squares  program  proved  extremely  useful  during  the  reanalysis  of 
earlier  data  as  it  permitted  repeated  reduction  of  the  same  experi¬ 
mental  data  with  incrementally  varied  input  parameters. 

The  FORTRAN  programs  together  with  a  glossary  of  terms  is  in¬ 
cluded  in  the  Appendix fA  thru  E.^ 

To  check  the  operation  of  the  FORTRAN  programs,  they  were  used 
to  perform  the  same  least  squares  and  solid  angle  correction  computa¬ 
tions  as  performed  earlier  by  Verser  and  Conway,  and  the  results  com¬ 
pared  with  their  earlier  results,  as  follows. 

1.  The  same  geometrical  parameters  as  used  by  Verser  and 
Conway  were  used  as  the  entering  arguments  for  the 
FORTRAN  solid  angle  correction  program.  Good  agree¬ 
ment  was  reached  with  Verser* s  hand-calculated  values 
and  Conway's  computed  values. 

2.  The  raw  data  (scaler  counts)  or  Verser  were  used  as 

the  entering  arguments  for  the  FORTRAN  least  squares 

program.  Good  agreement  was  reached  with  Verser1 s 

values  of  the  mean  values  of  W  (9.),  and  with  his 

exx  5 

coefficients  for  the  Legendre  series  representation 
of  W^g(9)  thus  indicating  proper  operation  of  the 
preliminary  and  subsequent  portions  of  this  program. 
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3.  The  results  of  analysis  of  Conway* s  data  with  the 
FORTRAN  program  did  not  agree  with  the  results  of 
analysis  which  he  reported.  Some  significant 
arithmetical  errors  were  uncovered  in  the  hand 
calculations  which  he  was  required  to  perform 
prior  to  entering  his  program.  When  these  errors 
were  corrected,  agreement  was  realized  in  the  re¬ 
sults  of  analysis  by  the  two  programs,  providing 
another  check  of  the  FORTRAN  least  squares  program. 

G.  Reanalysis  of  Earlier  Data 

Once  assured  that  the  new  FORTRAN  programs  were  operating  pro¬ 
perly  s  Verser's  data  and  a  portion  of  Conway1 s  data  were  reanalyzed, 
including  in  the  calculation  of  the  variances  the  effect  of  the 
electronic  instability  of  the  coincidence  unit. 

In  checking  over  the  data  collected  by  Conway,  it  was  noted  that 
they  were  collected  during  two  relatively  short  periods  of  time  spaced 

nearly  four  months  apart.  During  these  four  months  a  significant  change 

60 

in  the  strength  of  the  Co  source  (T  ~5.2  years)  occurred;  this 
had  not  been  corrected  for  in  his  data  analysis.  The  reanalysis  of 
his  experiment  was  done  using  only  his  later  set  of  data. 

It  became  apparent  that  the  scatter  of  Conway's  data  about  W^g(9) 
was  greater  than  the  standard  deviations  would  imply.  Consequently 9 
estimates  of  the  values  of  (T^Wex(@^)j  were  calculated  based  upon  the 
actual  residuals  between  the  experimental  points  and  Wjg(9)  for  his 
later  set  of  data.  This  calculation,  which  took  into  account  the  elec¬ 
tronic  variations  of  the  equipment,  resulted  in  values  for 


W 

ex 

W 


ex 


(V 


approximately  three  times  as  large  as  the  values  which  Conway  quoted. 

The  results  of  the  reanalysis  of  the  previous  data  are  tabulated 
below  and  are  presented  graphically  in  Figure  6a  and  6b.  It  is  evident 
that  the  anomaly  in  the  120-135°  region  still  exists  in  Verser's  results, 
though  it  is  not  as  significant  statistically  as  previously  reported  due 
to  the  larger  standard  deviations.  The  analysis  of  Conway's  data  now 
shows  no  evidence  of  a  dip. 
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1.30- 


Reanalysis  of  Data  of  Verser  and  Conway 
Fig.  6a,  b 
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Table  1 

Results  of  Reanalysis  of  Earlier  Data 


Legen 

a 

o 

dre  Coefficients 

a2 

a4 

Theoretical 

1.0000 

0.0996 

0.0084 

Verser 

1.0000  +  0.0042 

0.0960  +  0.0068 

0.0337  +  0.0087 

Conway 

1.000  +  0.011 

0.126  +0.020 

0.032  +  0.023 

An: 

Lsotropy 

Theoretical 

0.162 

Verser 

0.171  ±  0.026 

Conway 

0.220  +  0.079 

The  possibility  was  investigated  that  errors  could  have  been  intro¬ 
duced  into  Verser1 s  data  analysis  through  the  use  of  incorrect  parameters. 
In  investigating  this  possibility,  raw  data  gathered  by  both  Verser  and 
Conway  were  used  as  the  entering  arguments  for  the  FORTRAN  least  squares 
program,  and  these  data  were  reduced  while  systematically  varying  the  ex¬ 
perimental  parameters  to  see  what  effect  this  would  have  on  The 

results  of  this  analysis  are  as  follows: 


Table  2 

Effect  on  W^CG)  of  Varying  Input  Parameters 


Parameter  Varied 

7o  Variation  in 

Resulting  %  Variation 

Parameter 

Legendre  Coefficients 

Anisotropy 

a2 

a4 

r 

+  10% 

+8% 

+8% 

+8% 

Background 

+  6% 

0% 

0% 

0% 

Scaler  Resolving  Time 

+500% 

+3% 

+3% 

+5% 

It  can  be  seen  from  the  above  that  the  shape  of  W^s(Q)  is  relatively  in¬ 
sensitive  to  changes  in  background  and  changes  in  the  resolving  time 
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of  the  scalers.  However  the  shape  of  W..  (0)  varies  appreciably 

l  s 

when  T  varies.  The  reason  for  these  changes  in  W^s(9)  is  clear 

from  the  following.  If  the  numerical  value  used  for  T  is  greater 

than  the  actual  T  of  the  coincidence  unit  during  the  experiment, 

then  too  many  accidental  coincidences  are  subtracted  from  the  total 

coincidences  in  calculating  the  number  of  genuine  coincidences. 

Thus  the  apparent  value  of  C  is  smaller  than  was  actually  realized 

G 

experimentally.  Since  the  error  introduced  into  is  constant  at 

all  angles,  assuming  that  and  are  independent  of  9,  which  is 

approximately  true,  the  error  introduced  into  C  is  also  constant 

G 

at  all  angles.  Since  in  these  experiments,  the  ratio  of  the  acci¬ 
dental  coincidence  counting  rate  to  the  genuine  coincidence  counting 
rate  is  about  2:3  and  it  can  be  seen  from  the  expression 

W  (180)  \ 

rar  - 1  -  • 167  (3) 

Is 


that  W-  (90)  is  more  sensitive  than  W.  (180)  to  changes  in  the  genuine 
Is  Is 

coincidence  counting  rate,  it  is  apparent  that  the  calculated  anisotropy 
increases  when  the  coincidence  resolving  time  is  increased.  It  is 
stressed  that  the  shape  of  W^g(0)  is  distorted  by  using  a  value  of  T 
in  the  data  reduction  which  does  not  correspond  to  the  value  of  T 
actually  characteristic  of  the  coincidence  unit  during  the  experi¬ 
ment.  The  shape  of  W..  (9)  is  independent  of  T  as  long  as  the  true 

Is  0 

experimental  value  of  T  is  used  in  calculating  C  . 

A 

The  effect  upon  the  shape  of  W..  (9)  of  using  incorrect  para- 

x  s 

meters  for  the  geometry  of  the  experiment  was  investigated  also,  by 
progressively  varying  the  entering  arguments  for  the  FORTRAN  solid 
angle  correction  program.  The  results  of  this  analysis  are  given  in 
in  Table  3. 
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Table  3 

Effect  on  Solid  Angle  Corrections  of  Varying  Input  Parameters 


Parameter  Varied 

Amount  Varied 

Relative  Ch 
Absorption 

ange  in 
Coefficients 

Linear  absorption 
coefficient  of  the 
crystal 

+  107, 

0.00% 

0 . 00% 

Distance  from  source 
to  crystal 

+  8.5% 

0.00% 

0 . 00% 

Crystal  radius 

+  1.0% 

0 . 00% 

o 

O 

O 

Crystal  thickness 

+  1.0% 

0 . 00% 

0.00% 

It  can  be  seen  that  the  shape  of  W^s(@)  is  relatively  insensitive  to 
changes  in  these  geometrical  parameters  over  the  range  investigated. 

It  is  concluded  from  these  investigations  of  the  effects  upon  the 
least  squares  analysis  and  solid  angle  correction  computations  of  varying 
the  entering  arguments,  that  none  of  these  can  account  for  the  dip  re¬ 
ported  by  Verser, 
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H  Analysis  of  New  Data 


The  results  of  the  analysis  of  the  data  obtained  with  the  old 
slow  coincidence  unit  in  experiments  in  two  separate  quadrants  are 
indicated  graphically  in  Fig,  7.  When  the  data  from  each  of  these 
two  experiments  were  reduced  separately,  W^g(9)  agre^d  with  Wt^(9) 
within  the  statistical  error  of  the  experiment.  Since  the  angular 
correlation  function  is  symmetrical  about  180°,  it  was  possible  to 
combine  the  results  of  the  analysis  of  the  data  obtained  during  the 
experiment  in  the  90  -  180°  quadrant  with  those  obtained  during  the 
the  experiment  in  the  180  -  270°  quadrant.  The  purpose  of  combining 
the  results  in  this  way  was  to  benefit  from  the  decrease  in  (f  Wex(£K) 
resulting  from  the  analysis  of  a  greater  quantity  of  data,  thus  being 
able  to  compare  more  precisely  W^(9)  and  Wt^(9) 

Due  to  the  time  interval  between  the  two  experiments,  during 

which  the  source  decayed  significantly,  the  raw  data  from  the  two 

experiments  cannot  be  combined.  Consequently,  the  data  obtained  from 

each  experiment  were  analyzed  independently,  then  the  normalized 

values  of  the  weighted  mean  of  ^x(9^)  from  the  two  experiments  were 

combined  using  as  weighting  factors  the  reciprocals  of  their  variances. 

These  composite  values  of  W  (9.)  were  then  used  as  the  entering 

ex  i 

arguments  for  the  least  squares  subroutine  to  generate  a  new  (9) . 

Is 

The  results  of  the  new  experiments  are  tabulated  below. 


Table  4 

Results  of  Analysis  of  New  Data 


i - 

Experimental 

Theoretical 

Quadrant 

Composite 

180°-270° 

90-180° 

Anisotropy 

0.151  +  0.029 

0.165  +  0.034 

0.158  +  0.018 

0.162 

Legendre 

Coefficients 

a 

o 

1.0000  +  0.0057 

1.0000  +  0.0057 

1.0000  +  0.0040 

1.0000 

a2 

0.092  +  0.010 

0.1062  +  0.0094 

0.0987  +  0.0069 

0.0996 

a. 

0.010  +  0.012 

-0.005  +  0.012 

0.0036  +  0.0085 

0.0084 

4 

— 

i 
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Results  of  Analysis  of  New  Data 
Fig.  7. 
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I„  Comparison  of  Experimental  Results  with  Theory 

All  of  Verser*s  data,  and  the  portion  of  Conway* s  data  used  in 
the  reanalysis  of  his  data  (see  Section  G)  were  combined  with  the 
data  obtained  during  the  two  new  experiments  according  to  the  scheme 
discussed  in  the  previous  section.  That  is,  the  data  from  each  of 
the  four  experiments  were  analyzed  independently,  then  the  weighted 
means  of  the  normalized  mean  values  of  Wex(9..)  from  the  four  experi¬ 
ments  were  used  as  the  entering  arguments  for  the  least  squares  sub¬ 
routine.  The  results  of  this  combination  of  these  data  are  indicated 
graphically  in  fig*  8,  and  are  tabulated  in  Table  5»  No  significant 
deviation  from  Wt^(9)  is  evident. 

While  the  real  demonstration  of  the  proper  functioning  of  the 
experimental  equipment  lies  in  the  results  of  the  new  experiments  in 
two  quadrants,  it  is  of  interest  that  no  significant  deviation  from 
Wt^(9)  occurs  even  when  Verser's  and  Conway* s  results  are  combined 
with  these.  It  is  concluded  that  the  dip  which  they  reported  is 
strictly  a  statistical  fluke  and  that  there  is  no  actual  physical 
effect  involved. 


Table  5 


Results  of  Composite  Data  Analysis 
(Composite  Results  based  upon  4  experiments) 


- - , 

Experimental 

- f 

Theoretical 

_  _  ) 

Anisotropy 

t-  •  ~~  • 

0.169  +  0.018 

0.162 

Legendre 

Coefficients 

a 

1.0000  +  0.0028 

1.0000 

o 

a^ 

0.1000  +  0.0047 

0.0996 

2 

— 

a. 

0.0183  +  0.0058 

0.0084 

4 
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Results  of  Composite  Data  Analysis 


Fig.  8. 


26 


Ill 


DEVELOPMENT  OF  NEW  INSTRUMENTATION 

A.  New  Slow  Coincidence  Unit 

While  investigating  the  earlier  reported  anomaly  in  W^(@)  in 
the  105-135°  region,  work  proceeded  on  building  and  testing  a  new 
slow  coincidence  unit  incorporating  certain  portions  of  the  design 
reported  by  Garg  13^  .  This  design  features  a  variable  resolving 
time  which  is  independent  of  the  counting  rate.  The  slow  coinci¬ 
dence  unit  used  by  Verser  and  Conway  lacked  both  of  these  features 0 
The  coincidence  unit  which  was  constructed  differed  markedly 
from  Garg's  only  in  that  the  separate  channel  amplitude  discrimi¬ 
nators  of  Garg's  unit  were  deleted  since  differential  discrimina¬ 
tion  was  performed  earlier  (Fig.  9  ).  Figure  10  gives  the  com¬ 
plete  circuit  diagram  for  the  unit  which  was  constructed. 

The  new  slow  unit  was  tested  by  incorporating  it  into  the 
circuitry  used  for  the  experiment  in  the  90-180°  quadrant  in 
parallel  with  the  slow  coincidence  unit  which  Verser  and  Conway  had 
used.  The  correlation  function  measured  by  the  two  slow  coinci¬ 
dence  units  agreed  within  the  uncertainties  in  their  resolving  times. 
The  results  of  numerous  resolving  time  tests  of  the  new  slow  coinci¬ 
dence  unit  indicated  a  greater  long  time  drift  than  is  characteristic 
of  the  old  slow  unit.  This  increased  drift  of  T  necessitates  more 
frequent  determinations  of  its  magnitude  and  required  using  different 
values  for  T  during  the  angular  correlation  runs,  rather  than  using 
a  mean  value  of  f  based  upon  a  large  number  of  resolving  time  runs 
such  as  was  done  with  the  old  slow  coincidence  unite 

B.  Requirement  for  a  Fast  Coincidence  Unit 
I.  General 

A  coincidence  unit  of  significantly  shorter  resolving  time 

than  that  of  the  unit  used  by  Verser  and  Conway  was  required  to 

improve  the  statistical  error  in  the  experimental  estimate  of 

the  correlation  function,  W  (9.).  Better  statistics  were  re- 

ex'  i' 

quired  for  the  reinvestigation  of  the  earlier  reported  anomaly 
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Figure  9o 


Block  Diagram,  New  Slow  Coincidence  Unit 
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Figure.  10  a.  Circuit  Diagram  of  New  Slow  Coincidence  Unit 
Initial  Stages 
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Figure  10  b  .  Circuit  Diagram  of  New  Slow  Coincidence  Unit 
Final  Stages 
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2. 


in  the  value  of  W  (0.)  for  Ni^  in  the  105-135°  region,  and 

for  future  experiments  with  short-lived  nuclides  produced  in 

the  US  Naval  Postgraduate  School  Van  de  Graaff  generator. 

It  was  not  possible  to  reduce  the  resolving  time  of  the 

slow  coincidence  unit  below  a  value  T  .  which  is  determined 

mm 

by  the  width  of  the  window  of  the  differential  pulse  height 
analyzer  which  precedes  it.  This  time  corresponds  to  the  dif¬ 
ference  in  time  between  the  crossing  of  the  lower  voltage  dis¬ 
crimination  level  of  a  pulse  of  the  lowest  permissable  amplitude 
in  one  channel  and  the  time  of  crossing  of  the  lower  discrimi¬ 
nator  level  of  a  simultaneous  pulse  of  the  highest  permissable 
voltage  in  the  second  channel.  Because  of  this  irreducible 
minimum  value  for  T  y  and  since  a  resolving  time  larger  than 
7\n  gives  too  many  accidental  coincidences  for  an  acceptable 
variance  of  the  experimentally-determined  correlation  function, 
a  different  mode  of  operation  was  required.  The  scheme  used 
consisted  of  adding  a  fast  coincidence  unit  in  parallel  with 
the  slow  unit  and  mixing  the  output  pulses  from  the  two  units; 
this  scheme  reduced  the  accidental  coincidence  counting  rate 
because  of  the  decreased  T  of  the  fast  unit,  yet  retained  the 
amplitude*  discriminating  characteristics  of  the  differential 
pulse  height  analyzer  employed  with  the  slow  unit  (Figure  11  <>) 
The  fast  coincidence  unit  was  constructed  after  a  design 
by  Draper  and  Fleischer  (14)  and  was  adjusted  for  a  resolving 
time  of  approximately  20  nanoseconds,  thus  decreasing  C  by 

A 

about  a  factor  of  10  from  that  obtained  with  the  slow  units. 
Relation  between^  W^x(@^)j  and 


It  is  desired  to  minimize 

r  I 


(1) 


W  (0.) 
exv  i7 


Since  W  (0 . )  = 
exv  i 


N1N2 


and  the  fractional  deviations  of  N, 
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Figure  11 .  Block  Diagram  of  Experimental  Setup 
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are  very  small  compared  to  the  fractional  deviation  of 


<T 


“ex<6l>1  . 


(12) 


w  (0.) 

ex  1/ 


Effect  on  (T  (C  )  of  Decreasing  the  Coincidence  Resolving  Time 
G 

The  variance  of  C  is  calculated  from  equation  H  , 

G 


2  *  T 

(T  (CG)  =  -f 


<n. 


elec 

r 


(Jh 


•  2 
CA  + 


t 


It  is  evident  from  the  relation  C  =  2TN-,N~  (equation  6 1  that 

A  12|  I 

decreasing  the  resolving  time  of  the  coincidence  unit  7]  by  an 
order  of  magnitude  results  in  a  similar  decrease  in  C^.  This 
decrease  is  reflected  as  a  decrease  in  the  third  term  of  equation 
11  by  an  order  of  magnitude  since  t  is  constant  in  these  experi¬ 
ments  o  Assuming  no  change  in  the  fractional  resolving  time  fluctu 
“eieM 

the  second  term  of  equation  11  would  decrease 
The  genuine  coincidence  counting  rate 


ation,  \  j3 

by  a  factor  of  about  100, 


is  independent  of  the  coincidence  unit  resolving  time,  as  long 
as  the  latter  exceeds  the  average  interval  between  the  emissions 
of  the  two  quanta;  this  is  the  case  for  all  three  coincidence 


units  discussed  in  this  report . 


Thus  a  decrease  in  C  also  re~ 

A 

This  decrease  in 


suits  in  a  decrease  in  C  since  C  =  C  +  C  , 

T  £  T  A  G 

C  results  in  a  decrease  in  _T.  Therefore  all  three  terms  of 

T  t 

equation  11  decrease  when  T  decreases. 

Relative  Effect  of  Coincidence  Resolving  Time  Drift 

2  • 

The  previous  section  indicated  how  (j~(C  )  is  decreased  by 

G 

decreasing  T  .  Then  from  the  relation 


<r 


tHe*<y] 


<T<cr> 


equation  12 , 
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it  is  evident, 


since  is 


constant, 


that 


<r 


Kx<°l>] 

UCx<°l> 


has  been  decreased  appreciably. 


The  assumption  that 


,<wr  >J 


did  not  vary  was  made  in  the  last  section  to  simplify  the  dis¬ 
cussion.  It  is  evident,  however,  that  due  to  the  large  decreases 

in  the  values  of  all  other  terms  in  equation  11 ,  including  the 
*  2 

term  which  multiplies 
latter  can  be  tolerated  while  still  decreasing  very  significantly 
the  fractional  variance  of  Wex(CL) *  Thus  the  electronic  stability 
of  the  fast  coincidence  unit  is  not  as  critical  as  it  would  be 
with  a  unit  of  larger  coincidence  resolving  time. 


^lec^  ^  ' 

r 


,  larger  values  of  the 


Co  Operation  of  the  Fast  Coincidence  Unit 

The  basic  theory  of  operation  of  the  fast  unit  and  the  circuitry 
for  the  Draper  Fleischer  unit  are  given  in  detail  elsewhere  |^14]  and 
will  not  be  repeated  heref  Either  single  channel  becomes  operative 
when  the  input  pulse  reaches  68%  of  its  maximum  amplitude,  regardless 
of  the  value  of  that  amplitude  within  the  range  20  100  volts.  The  cir¬ 
cuit  is  so  designed  that  resolving  times  on  the  order  of  100  times 
shorter  than  the  rise  time  of  the  input  pulses  can  be  achieved.  This 
characteristic  is  required  since  the  rise  time  of  the  pulses  from  the 
sodium  iodide  scintillation  crystals  is  on  the  order  of  200-300  nano¬ 
seconds  o 

Only  the  coincidences  of  1.17  MEV  and  1„33  MEV  gamma  rays  are  of 
interest  in  this  experiment  with  Ni^.  Amplitude  discrimination  which 
was  performed  by  the  differential  pulse  height  analyzers  which  preceded 
the  slow  coincidence  unit  selected  the  pulses  corresponding  to  the 
photopeaks  of  the  two  gamma  rays.  All  coincidences,  regardless  of  the 
pulse  height,  were  detected  by  the  fast  unit  since  no  amplitude  dis¬ 
crimination  preceded  it.  The  output  pulses  of  the  slow  coincidence 
unit  and  those  of  the  fast  unit  were  passed  through  a  mixer  circuit 
which  effectively  functioned  as  another  coincidence  unite  An  output 


*  A  block  diagram  of  the  fast  unit  is  provided  in  Fig.  12. 
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Block  Diagram,  Fast  Coincidence  Unit 
Figure  12. 


*  After  Farley^  [15] 

**  After  Bell  &  Petch  &  Graham  ,  M 


pulse  was  obtained  from  the  mixer  only  if  there  was,  in  a  time  interval 
of  about  35  microseconds,  a  20  nanosecond  coincidence  between  pulses 
of  any  height  a$  ascertained  by  the  fast  unit  and  a  200-nanosecond 
coincidence  between  pulses  corresponding  to  the  photopeaks  of  interest 
as  determined  by  the  slow  unit.  Thus  by  combining  these  two  units 
in  this  manner,  an  energy-selective  counting  system  with  a  20-nano- 
second  coincidence  resolving  time  which  was  independent  of  the  counting 
rate  was  obtained. 

Certain  modifications  were  necessary  to  the  fast  unit  and  are 
given  in  Figure  13,  The  last  two  stages  of  the  Draper-F leischer 
unit  could  not  be  made  to  function  properly.  This  portion  of  the 
circuitry  was  therefore  taken  out,  and  a  cathode  follower  was  sub¬ 
stituted  . 

The  mixer  circuit  which  contained  circuitry  for  amplitude  dis¬ 
criminating,  shaping,  amplifying,  and  mixing  the  fast  coincidence 
pulse  with  the  slow  coincidence  pulse  was  built  up  and  tested  stage- 
by-stage.  The  discriminator  of  the  mixer  unit  is  a  uni-stable  multi¬ 
vibrator  with  Zener  diode  stabilization.  The  inverter  stage  in  the 
fast  coincidence  channel  of  the  mixer  results  in  a  negative  pulse 
of  approximately  35-microsecond  duration,  the  duration  varying  with 
the  discriminator  setting.  This  pulse  is  then  mixed  with  a  0.7  micro¬ 
second  negative  pulse  from  the  slow  coincidence  unit.  Discrimination 
of  the  mixed  output  is  performed  by  the  scaler  counting  these  fast- 
slow  coincidences.  The  design  of  the  mixer  circuit  is  given  in  Fig.  14. 

Do  Testing  and  Alignment  of  the  Fast  Coincidence  Unit 

The  alignment  of  the  fast  coincidence  unit  consisted  of  two 
major  portions-achieving  linearity  in  each  channel,  and  assuring  the 
detection  of  true  coincidences  by  the  coincidence  stage. 

Draper  and  Fleischer  ^lJ  describe  a  means  of  testing  for 
linearity;  i.e.,  ensuring  that  the  fast  discriminator  fires  when  the 
input  reaches  the  same  predetermined  fraction  of  its  full  amplitude, 
over  a  wide  range  of  amplitudes.  They  prescribe  a  mercury-relay  pulser 
delivering  a  five-volt  negative  step  with  a  rise  time  of  seven  nano¬ 
seconds  which  is  split,  one  portion  triggering  a  10-mc  ringing  cir¬ 
cuit  which  provides  the  time  reference,  the  other  going  through  an 
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Since  pulses  from  the  fast  discriminator  (V9,  V10  of  the  Draper-Fleischer  circuit)  are  not  square, 
it  is  necessary  to  amplitude  discriminate  the  output  of  the  fast  coincidence  unit.  This  is  accomplished 
by  amplifying  the  mixed  output  of  the  two  time  compensation  units  in  and  and  then  amplitude  dis¬ 
criminating  in  the  circuit  comprised  of  V3  and  V^.  The  output  of  V4  then  is  inverted  and  mixed  with  the 
output  of  the  slow  coincidence  unit  (Figure  14  ). 


Figure  13.  Fast  Coincidence  Pulse  Amplifying  and 
Discriminating  Circuitry 
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Figure  14. 


Circuit  Diagram  of  Mixer  Circuit 


attenuator,  linear  pulse  amplifier,  and  the  coincidence  unit 
(Figure  15 -a) . 

Several  modifications  of  this  test  scheme  were  tried,  due  to 
the  nonavailability  of  a  pulse  generator  which  could  deliver  a  pulse 
of  sufficiently  large  amplitude  and  short  rise  time  to  trigger  the 
ringing  circuit  satisfactorily,  and  due  to  the  lack  of  a  compensated 
attenuator.  These  modifications  met  with  varying  degrees  of  success. 
When  a  double  pulse  generator  was  used  as  a  signal  source,  and  the 
amplitudes  of  the  inputs  to  the  separate  channels  of  the  fast  unit 
were  varied  by  adjusting  the  gains  of  the  linear  amplifiers,  dis¬ 
tortion  of  the  wave  shapes  resulted  when  the  coarse  gain  settings 
were  varied.  This  made  it  impractical  to  vary  the  amplitudes  of 
the  coincidence  unit  inputs  by  varying  the  gains  of  the  linear 
amplifiers  since  the  time  reference  for  determining  the  delay  within 
the  coincidence  unit  itself  was  thus  lost.  Still  greater  wave  shape 
distortion  resulted  from  varying  the  gain  of  the  double  pulser. 

The  procedure  which  has  finally  emerged  as  workable  is  outlined 
in  Figure  15-b  and  consists  of  triggering  the  oscilloscope  with 
the  output  of  the  fast  unit  discriminators  and  viewing  the  relative 
delay  of  nuclear  pulses  of  varied  amplitudes.  The  problem  of  wave 
shape  distortion  when  the  coarse  gain  settings  of  the  linear  ampli¬ 
fiers  were  adjusted  has  been  avoided  by  limiting  the  adjustment  of 
delay-equivalence  to  the  10-volt  range  of  amplitudes  used  experi¬ 
mentally.  Since  a  10-volt  window  is  used  in  the  differential  pulse 
height  analyzer  experimentally,  an  equivalent  time  delay  over  this 
range  is  all  that  is  required.  The  inputs  to  the  fast  unit  are  varied 
over  this  10-volt  range  by  adjusting  only  the  fine  gain  controls, 
avoiding  the  distortion  of  wave  shape  which  accompanies  changes  of  the 
coarse  gain  settings.  The  adjustment  of  the  fine  gain  controls  can 
be  avoided  by  providing  a  nuclear  source  or  a  combination  of  sources 
with  different  decay  energies  which  results  in  different  amplitude 
signals  being  received  by  the  fast  unit  simultaneously.  The  relative 
delay  introduced  into  these  pulses  of  different  amplitudes  is  then 
observable  directly. 
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Fig,  15  o  Fast  Coincidence  Unit  Test  Scheme 
(a)  (As  described  by  Draper  and  Fleischer) 


(b)  As  employed 
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An  alternate  means  of  accomplishing  this  alignment  consists  of 
using  test  pulses  from  a  double  pulser,  triggering  the  oscilloscope 
with  the  trigger  output  pulse  of  the  double  pulser  and  viewing  the 
output  of  the  fast  discriminator  of  the  fast  coincidence  unit  in 
each  channel.  Again,  the  input  amplitude  is  varied  by  changing  the 
settings  of  the  fine  gain  controls  on  the  linear  amplifiers. 

The  procedures  already  described  ensure  only  that  the  delay 
which  a  pulse  experiences  as  it  passes  through  the  time  compensating 
circuit  of  one  channel  is  independent  of  amplitude „  The  other  major 
alignment  requirement,  which  cannot  be  met  until  each  channel  has 
been  aligned  separately,  is  the  assurance  that  the  true  coincidences 
are  being  detected  by  the  unit.  This  requires  identical  total  delay 
from  the  time  of  emission  from  the  source  through  all  the  circuitry 
associated  with  each  of  the  two  channels  and  through  the  fast  coinci¬ 
dence  unit  itself  up  to  the  coincidence  stage. 

To  determine  the  ability  of  the  fast  unit  to  detect  true  coinci¬ 
dences  before  it  was  completely  aligned,  it  was  used  in  parallel 
with  a  slow  unit  in  the  angular  correlation  experiment  in  the  90-180° 
quadrant.  While  it  generally  tracked  the  slow  unit,  it  did  so  errati¬ 
cally. 

Variable  delay  lines  have  been  installed  to  facilitate  the  align¬ 
ment  for  equivalent  delay  in  the  two  signal  channels  to  the  coincidence 

stage  of  the  fast  unito  The  alignment  is  accomplished  using  the  experi- 
60 

mental  Co  source  with  the  experimental  geometry  and  normalizing  the 
total  coincidence  counting  rate  measured  by  the  fast-slow  coincidence 
detecting  scheme  to  the  genuine  coincidence  counting  rate  determined 
by  the  slow  coincidence  unit.  After  normalizing  the  total  coincidence 
counting  rate  of  the  fast-slow  combination  in  this  manner,  and  making 
a  graph  of  this  versus  the  setting  of  the  variable  delay  line  in  one 
of  the  time  compensation  channels  of  the  fast  unit,  a  plateau  is  ob 
tained  the  width  of  which  corresponds  to  the  resolving  time  of  the 
fast  coincidence  unit.  The  delay  line  is  then  adjusted  to  the  center 
of  this  plateau  as  determined  from  the  graph. 
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Thus  procedures  have  been  developed  for  aligning  the  time  com 
pensation  circuit  of  each  channel  for  linearity  by  means  of  either 
a  nuclear  source  or  a  double  pulser,  and  for  aligning  the  two  chan¬ 
nels  of  the  fast  unit  for  equivalent  delay  by  means  of  a  nuclear 
source . 

E  Equipment  Status 

The  experimental  equipment  described  herein  including  the  new 
fast'  slow  coincidence  detecting  scheme  is  functioning  well  enough 
to  be  used  in  the  experiments  for  which  it  was  assembled;  namely, 
the  angular  correlation  of  short'  lived  nuclides  activated  in  the 
Van  de  Graaff  generator.  Due  to  limited  unavoidable  drifts  in  the 
electronic  equipment,  frequent  checks  must  be  made  of  operating 
characteristics,  and  appropriate  adjustments  made  as  necessary  to 
ensure  the  closest  possible  control  over  the  experimental  environ¬ 
ment.  The  checkout  procedures  which  have  proved  extremely  useful 
in  ensuring  that  easily'  overlooked  actions  are  not  overlooked  and 
that  checks  are  made  in  the  same  sequence  for  the  greatest  possible 
reproducibility  of  procedures,  should  be  continued 

Figure  16  illustrates  the  layout  of  the  major  electronic  portion 
of  the  experimental  equipment,  all  rack-mounted  and  mobile  to  facili¬ 
tate  its  use  with  or  without  the  Van  de  Graaff  generator. 
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Fig.  16  Electronic  System 


FORTRAN  PROGRAM 


Appendix  Part  A 


Least  Squares  Analysis 


44 


PROGRAM  ANGCORR 

DIMENSION  AN  1 ( 20 »  20 ) ,  AN2(20,20),  CT(20,20),  83(20), 

1  K ( 20 ) ,  AND  I P ( 20, 20 ) ,  AND2 P ( 20 , 20 ) ,  CTD(20,20),  CAD(  20,20), 

2  W  (  20 , 20 )  ,  SW  (20,20)  ,  BK20),  82(20),  SBB(20),  D(20),  ANIS0(2) 

3  ,TAUC( 20,20) , STAUC ( 20, 20 ) , TAUD 1 ( 20 , 20 ) , TAUD2 ( 20 , 20  ) ,  BACK  1(20, 20 

4  ) , BACK2( 20, 20), T( 20, 20) ,  CGD(20,20),  JRUN(20,20) 

COMMON  Bl,  B2 ,  B3,  L,  ALPHA1,  ALPHA?.,  ALPHA3,  SBB 
NN  =  1 

4  FORMAK  3F  10.6) 

6  FORMAT ( 4 1 HOGAMMA-GAMMA  ANGULAR  CORRELATION  ANALYSIS/) 

7  FORMAT ( 38H0NUMBER  OF  SETS  OF  DATA  TO  BE  ANALYZED/) 

WRITE  OUTPUT  TAPE  4,  6 

398  FORMAT ( 23H023  APRIL  62  ANG  CORR  ) 

WRITE  OUTPUT  TAPE  4,  398 
WRITE  OUTPUT  TAPE  4,  7 

8  FORMAT  (12) 

9  FORMAT  (110) 

10  FORMAT  (F12.ll,  3F12.10,  F7.4,  F8.4,  F9.2) 

11  FORMAT ( F6. 2 ,  110) 

12  FORMAT  (F9.1,  F20.1,  F17.1,  134) 

13  FORMAT( F 10.2,  13/) 

14  READ  INPUT  TAPE  3,  8,  MM 


q 

500  FORMAT ( 13//) 

WRITE  OUTPUT  TAPE  4,  500,  MM 


O' 


501  FORMAT! 13) 

15  READ  INPUT  TAPE  3,  501,  L 

READ  INPUT  TAPE  3,  1 0 , CT AUC , CST AUC , CTAUD 1 , CT AUD2 , C8ACK 1 , C8 ACK2 , AT 

READ  INPUT  TAPE  3,  11,  ( 83 ! J )  ,  K(J),  J  =  1,1) 

REAO  INPUT  TAPE  3,  4,  ALPHA1,  ALPHA2 ,  ALPHA3 

16  FORMAT! 16HOOATA  SET  NUMBER/) 

WRITE  OUTPUT  TAPE  4,16 
WRITE  OUTPUT  TAPE  4,  500,  NN 
REAO  OUT  DATA  INPUT  FOR  CHECK 

516  FORMAT! 23H0CHECK  INPUT  PARAMETERS/) 

WRITE  OUTPUT  TAPE  4,  516 

517  FORMAT 1 69H0  TAUC  STAUC  TAU01  TAU02  BACK  1 

1  BACK2  T/) 

WRITE  OUTPUT  TAPE  4,  517 

518  FORMAT  IF13.11,  3F12.10,  F7.4,  F8.4,  F9.2//) 

519  WRITE  OUTPUT  TAPE  4 , 5 1 8, C TAUC , CSTAUC , CT AUD 1 , CT AUD2 , C BACK  1 , C8 ACK2 , 
1  AT 


520  FORMAT ( 7 1HOTHEORETICAL  LEGENDRE  POLYNOMIAL  COEFFICIENTS  CORRECTED 


IFOR  SOLID  ANGLE// , 3F 12.6/ ) 

WRITE  OUTPUT  TAPE  4,  520,  ALPHA  1 ,  ALPHA2 ,  ALPHA3 
521  FORMAT (29H0NUMBER  OF  RUNS  AT  EACH  ANGLE/) 

WRITE  OUTPUT  TAPE  4,  521 

WRITE  OUTPUT  TAPE  4,  13,  !B3(J),  K(J),  J  =  1,  L) 

DO  18  J  =  1,  L 
II  =  1 

17  READ  INPUT  TAPE  3,  12,  A N  1  C  I  I , J ) »  AN2  C 1 1 , J ) ,  CT (  I  I , J ) , JRUN ( I  I , J ) 
READ  INPUT  TAPE  3,  10,  TAUCIII.J),  STAUC(II,J),  TAUD 1 (I  I , J ) , TAUD2I 
1  II, J),  BACK  1  ( 1 1 , J ) ,  8AC  K2 ( 1 1 , J ) ,  T(II,J) 

399  I F { T AUC ( II ,J 1-0.0)400,400,401 

400  TAUC 1 1 1 , J )  =  CTAUC 

401  IFISTAUCl I  I, J)-0.0)  402,  402,403 

402  S T AUC ( I  I , J )  =  CSTAUC 

403  IFITAUD1 (  I  I , J  )-0 . 0 ) 404 ,4 04 , 4  05 

404  TAUD I ( 1 1 , J )  =  CTAUDI 

405  I  F  C TAUD2I  I  I , J  )-0 . 0 ) 406 ,4 06 ,4  07 

406  TAUD2 ( 1 1 , J )  =  CTAUD2 

407  IF (BACK  I  ( I  I , J )-0 . 0 ) 408 , 4 08 ,4 09 


408  BACK  1 (  II, J)  =  CBACK1 


409  I F ( BACK2 ( II, J)-0. 0)410, 410, 411 

410  BACK2 ( I  I , J )  =  CBACK2 

411  IF(T(II,J)-0. 0)412, 412. 413 

412  T( II , J 1  =  AT 

413  II  =  II+1 

IF  (II-K(Jl)  17,  17,  18 

18  CONTINUE 

19  FORMAT (56H0C HECK  DATA  INPUT  OF  Nl,  N2,  AND  CT  AND  INPUT  PARAMETERS 


1/) 

WRITE  OUTPUT  TAPE 

4,  19 

520  FORMAT ( 9H0ANGLE- 

,  F6.2 / ) 

JJ=  1 

529  WRITE  OUTPUT  TAPE 

4,  528, 

B3( JJ) 

20  FORMAT ( 1 15H0  Nl 

N2 

CT 

TAUC 

STAU 

4> 

00 

1C  TAUD1 

TAUD2 

BACK1 

BACK2 

T 

RUN  /) 

21  FORMAT ( F 1 0. 1 , F 1 2. 1 , F 1 1 . 1 , F 1 3 . 1 1 , 3F 1 2. 1 0 , F7 . 4 , F8. 4 , F9 .2 , 19 ) 

WRITE  OUTPUT  TAPE  4,  20 
II  =  1 

530  WRITE  OUTPUT  TAPE  4,  21,  AN1(II,JJ),  AN2(II,JJ),  CT  (  1 1  ,JJ  )  ,  TAUC  (  I 
1  I  ,  JJ  ) ,  STAUCI  II  ,  JJ  )  ,  TAUD1  (II,  J  J  )  ,  TAUD2  (  1 1  ,  J  J  )  ,  BACK  1  (  I  I,  JJ);,  BACK2(  II 

•2  ,  JJ  )  ,T(  I  I  ,  JJ)  ,  JRUNUI,  JJ) 

I F ( I  I-K ( J J ) )  531,532,532 

531  11=11+1 
GO  TO  530 


532  FORMAT ( SHO  /) 

WRITE  OUTPUT  TAPE  4,  532 
IF(  JJ-U533,  534,534 

533  JJ=JJ+1 
GO  TO  529 

COMPUTE  NDOTI,  ND0T2 ,  CTDOT,  W  AND  SIGMA  W  FOR  EACH  RUN 

534  FORMAT  ( 20H0RESULTS  OF  ANALYSIS/) 

WRITE  OUTPUT  TAPE  4,  534 

JJ  =  I 
25  1 1 =  I 

30  AND1P(  I  I »  J  J  )  =  AN1 ( 1 1 , JJ ) /T< 1 1 , JJ )*( 1.0+TAUD1 ( 1 1 , JJ ) *AN1 ( 1 1 , JJ )/T 
I  (  I  I  ,  JJ  )  ) 


AN02P ( II, JJ)  =  AN2(II,JJ)/T(II,JJ)*( 1 ♦ 0  +  TAU02 ( I  I , J J ) *AN2 ( I  I , J J ) / T 
1  (II, JJ)) 

CTO ( 1 1 , J J )  =  CT(II,JJ)/T(II,JJ) 

CAD ( 1 1 , J J )  =  2.0*TAUC( II  , JJ)*AND1P( I  I, JJ)*AND2P( II, JJ) 

CGD  (  1 1 , J J  )  =  CTD ( I  I , J J )  -  CAD ( I  I , J J ) 

W(II,JJ)  =  ( CTD( I  I , JJ )-CAD(  I  I  ,  JJ)  )/ ( ( AND  1  P  C I  I , JJ J-BACK1 ( 1 1 , JJ ) ) 

1  * ( AND2P ( II»JJ)-BACK2(II »JJ) )) 

SWIII, JJ)  =  W( II , JJ)*SQRTF ( ( CTD ( 1 1 , J J ) /T ( 1 1 , JJ )  + ( CAO ( 1 1 , JJ ) /T ( 1 1 , J 
1  J))  +  ( STAUC( 1 1 , JJ)*CAD ( I  I , JJ )/TAUC( I  I , J J ) ) »*2 ) )/(CTD( I  I, JJ )-CAD 


2  (  II , JJ  )  ) 

IF  ( I I-K  C JJ)  )  35,  40,  40 
35  II  =  II  +  l 
GO  TO  30 

40  IF  (JJ-L)  45,  50,  50 
45  JJ  =  JJ+1 
GO  TO  25 

C  COMPUTE  MEAN  VALUES  OF  W,  SIGMA  W,  AND  RMS  SIGMA  FOR  EAIH  ANGLE 

50  JJ  =  1 
55  SUM  1  =  0.0 
SUM2  =  0.0 

II  =  1 

60  SUM  1  =  SUMJ+W(II,JJ)/(SW(II,JJ)**2) 

SUM2  =  SUM 2+ l.0/(SW(II,JJ)**2) 

IF  (I  I — K ( JJ )  )  65,  70,  70 
65  II  =  II  +  l 
GO  TO  60 

70  BHJJ)  =  SUM  I / SUM2 
B2(JJ)  =  SQRTF ( 1 . 0/SUM2) 

SUM3  =  0.0 

II  =  1 

71  SUM3  =  SUM3+ (W( I  I , JJ )-B1  ( JJ )  )**2 
IF  (II-K(JJ))  72,  73,  73 

72  II  =  II+l 


GO  TO  71 

73  D(JJ)  =  K  (  J J ) 

SBB(JJ)  =  SQRTF ( SUM3/( (D (JJ)-1»0)*D(JJ) ) ) 
IF  (JJ-L)  75,  80,  80 


75  JJ  =  JJ+1 
GO  TO  55 

C  READ  OUT  RESULTS  SO  FAR 

80  FORMAT  ( 11 8HONOOT 1 PRI ME  ND0T2PRI ME  CTOTALDOT 

1  COOT  CGENDOT  W  SIGMAW  RUN 

2) 

82  FORMAT  { F 1 0. 4 , F20 . 4, 2F 13 . 4 , F 1 4. 4 , 2E20. 5 , 1 5 ) 

WRITE  OUTPUT  TAPE  4,  80 
DO  87  J=  1,1 

84  FORMAT ( 9H0ANGLE-  ,  F6.2/1 

•  WRITE  OUTPUT  TAPE  4,  84,  B3(J) 

II  =  1 

85  WRITE  OUTPUT  TAPE  4,  82,  AND1P(II,J),  AND2P(II,J),  CTD(II,J), 

1  CAD ( 1  I , J  1 ,  CGD ( I  I , J ) ,  W(II,J),  SW(II,J),  JRUN (  1 1 , J  ) 

II  =  II+l 

IF  ( I I-K { J  11  85,  85,  87 
87  WRITE  OUTPUT  TAPE  4,  532 


CAC 

// 


90  FORMAT  ( 8 1  HO  ANGLE  WBAR  SI( 

IMA  WBAR  WBARB  /) 

WRITE  OUTPUT  TAPE  4,  90 
95  FORMAT  (F20.2,  3E20.5/J 

100  WRITE  OUTPUT  TAPE  4,95,(B3(J),  BHJ),  B2(J),  SBB(J),  J=1,L) 

C  START  LEAST  SQUARES  ANALYSIS 

105  CALL  LEASTSQ 
115  IF  (NN  -  MM)  120,  125,  125 

120  NN  =  NN+1 
GO  TO  15 

125  FORMAT! 15H0END  OF  PROBLEM) 

WRITE  OUTPUT  TAPE  4,  125 

END  FILE  4 

PAUSE 

END 

SUBROUTINE  LEASTSQ 
JK  =  1 

C  SET  UP  ARRAYS  AND  COMPUTE  LEGENDRE  POLYNOMIAL  MATRIX 

201  DIMENSION  P0LY{20,3),  OMEGA ( 20,20 ) ,  PART1(3,20),  CEE!3,3), 

1  CEE  I NV ( 3,3),  XI  13) ,  C0EFT13),  SC0EFT!3),  ARGP(20),  B 1 ( 20 ) , 


2  B2 ( 20  ) ,  B3! 20 ) ,  SUM713),  ANIS0(2),  SBB!20) 


COMMON  81,  82,  B3,  L,  ALPHA!,  ALPHA2,  ALPHA3,  5BB 

WFUNCTF( A,B,C,D)  =  A+ ( B» 0 . 25 ) » ( 3 . 0»C0SF { ( 2. 0*D ) /57 . 2957795 )  *  1 . 0 )  + 
1 (C*0.0 15625) *(35. 0*C0SF( (4.0*01/57.2957795)  +  20 .0*C0SF (( 2.0*0 ) / 

2  ' 57.2957795)  +  9.0) 

DO  203  1=1, L 

ARGP(I)  =  COSF( (B3(I )/57. 2957795) ) 

POLY ( I ,  I  )  =  1.0 

POLY (1,2)  =  0.5*(3.0*(ARGP( I ) **2 )-1 .0) 

203  POLY (1,3)  =  0. 125*135. 0* ( ARGP ( I ) **4 ) -30 . 0*(ARGP(  I)**2)+3.0) 
IF(JK-1 >206,206,204 

204  DO  205  I  =  1  »L 
POLY (1,1)  =  1.0 

POLY (1,2)  =  ARGP ( I ) *»2 

205  POLY (1,3)  =  ARGP ( I ) *»4 

C  COMPUTE  OMEGA  MATRIX,  PARTI  MATRIX  AND  CEE  MATRIX 

206  DO  210  1=1, L 
00  210  J= 1 »L 

210  OMEGA( I , J )  =  0.0 
DO  215  I  =  1,  L 

215  OMEGA (1,1)  =  ( (  1.0)/B2(I  ) )**2 
DO  225  1=1,3 
DO  225  J= 1 , L 
SUMS  =0.0 
00  220  K=1,L 


i 


220  SUM5  =  POL  Y I K , I  ) » OMEGA ( K  »  J )  +  SUM5 
225  PART  1  I  I ,  J  )  =  SUM5 
DO  235  1=1,3 
DO  235  J=1 ,3 
SUM6  =  0.0 
DO  230  K=  1  ,  L 

230  SUM6  =  PARTI  I  I ,K)*POLY(K , J)  +  SUM6 
235  CEE ( I ,  J  )  =  SUM 6 
INVERT  CEE 

OETCEE  =  C EE {  1  i 1  )  » (CEE ( 2 , 2 ) *CEE ( 3, 3 )-CEE ( 3 , 2 ) "CEE ( 2, 3 ) )-CEE(2, 1 ) 
1 CE  E (  1,2)»CEE(3,3)-CEE(3,2)*CEE( 1,3) )+CEE(3,l ) *  I  CEE  I  1 , 2 > *CEE( 2 , 3 > 
2C£E(2,2)*CEE( 1,3) ) 


CEEINVI  1»  1  )=(CEE(2,2)*CEE(3,3)-CEE(3,2)*CEE(2»3) ) /DETCEE 
CEEINVI  1»2)=(CEE(3»2)*CEE(  1,3) -CEE  I 1 , 2 ) *CEE ( 3 , 3 ) ) /OETCEE 
CEEINVI  1,3)  =  (CEE( 1 , 2 ) *CE E ( 2 , 3 >-CEE I 2 , 2 ) *CEE ( 1,3) J /DETCEE 
CEEINVI 2, 1 )  =  I  CEE  I  3 , 1 ) *CE  E I  2 , 3 ) -CEE (2,1 )*CEE(3,3) ) /DETCEE 
CEEINV(2,2)=(CEE( 1 , 1 ) *CE E ( 3 , 3 )-C EE ( 3 , 1) *CEE ( 1 , 3 ) ) /DETCEE 
CEEINV(2»3)=(CEE(2»1)*CEE( Ij3)-CEE( 1,1 )  »CEE I  2, 3 ) ) /DETCEE 
CEEINVI 3, 1 ) = I  CEE  I  2 , I ) *CE E I  3 , 2 )-CEE I  3 , 1 )*CEE(2,2) ) /DETCEE 
CEEINVI 3,2)  =  (CEE  I  3, 1) *CEE( 1 ,2)-CEE( 1 , 1 ) *C EE  I  3 , 2 ) ) /DETCEE 
CEEINVI 3,3)= (CEE  I  1 , 1 )*CEE (2,2)-CEE(2, 1 )*CEE(  1 ,2)  )/ DETCEE 


241  FORMAT  (3F20.12/) 

DO  246  1=1,3 
SUM83  =  0.0 
DO  245  K  =  l,L 

245  SUM83  =  PART  1  ( I , K ) *8 1 ( K)  +  SUM83 

246  XI(I)  =  SUM83 

COMPUTE  COEFFICIENTS  AND  THEIR  SIGMAS,  NORMALIZING  TO  FIRST 

COEFFICIENT  AND  COMPUTE  ANISOTROPY 

DO  248  1=1,3 

SUM81  =  0.0 

DO  247  K  =  1 ,  3 

247  SUMS  1  =  CEEINVl I ,K)*XI (K )  +  SUM81 

248  COEFTl I )  =  SUM81 
DIV  =  COEFTl 1 ) 

DO  249  1=1,3 

249  COEFTl I  )  =  COEFTIII/DIV 
DO  260  1=1,3 

260  SCOEFTl  I  )  =  SQRTF (CEE  I NV 1 1 , 1 ) ) /DI V 
AN  I  SO  (  2  )  =  C0EFT12)  +  COEFT  (3). 

265  IFIJK-2)  274,  340,  340 

274  AN  I  SO l  1  )  =  l WFUNCTF ( COEFT ( 1 ) ,  COEFT 1 2 ) ,  COEFT ( 3 ) ,  180.0)  - 

1  WFUNCTFICOEFTI  1  ) ,  COEFT ( 2 )  ,  COEFT ( 3 ) ,  090.0 )) /WFUNCTF ( COEFT (  1 ) , 

2  COEFT ( 2 ) ,  COEFT ( 3  ) ,  090.0)  ' 


275  FORMAT (71H0LEAST  SQUARE  FIT  COEFFICIENTS  OF  LEGENDRE  POLWOMIALS  A 
1  NO  THEIR  SIGMAS/) 

WRITE  OUTPUT  TAPE  4,  275 

280  WRITE  OUTPUT  TAPE  4,  342,  (COEFT(I),  SCOEFT(I),  I  =  1,3) 

292  FORMAT ( F9. 6/ ) 


293 


294 


295 


296 


300 


VEESQ1  =  0.0 
DO  293  J  =  1,L 

VEESQ1  =  (B1(J)/B2(J) )**2  +  VEESQ1 
VEESQ2  =  0.0 
DO  294  K  *  1,3 

VEESQ2  =  COEFT(K)*XI (K)*DI V  +  VEESQ2 
E  =  L 

EPSSQ  =  (VEESQ1  -  VEESQ2)/(E  -  3.0) 

FORMAT! 18H0EPSIL0N  SQUARE  =  ,  E20.6/) 

WRITE  OUTPUT  TAPE  4,  295,  EPSSQ 
IF! B3( 3 )- 1 80.0 )  296,  296,  320 

FORMAT! 45H0EXPERIMENTAL  FUNCTION  THEORETICAL 
COMPUTE  POINTS  ON  THEORETICAL  AND  EXPERIMENTAL 
WRITE  OUTPUT  TAPE  4,  296 
SUM9  =  90.0 

WFUNCT  =  WFUNCTF ( COEFT ! I ) ,  COEFT { 2 ) ,  COEFT { 3 ) , 


FUNCTION/) 

CURVES 


SUM9  )  /WFCNCTF 


315 

316 

1316 


317 

318 

320 


1 ( COEFT {  1 ) ,  COEFT {  2  )  ,  COE 
VFUNCT  =  WFUNCTF ( ALPHA  1, 
1  ALPHA2,  ALPHA3 ,  090.0) 
WRITE  OUTPUT  TAPE  4,  324 
SUM9  =  SUM 9  +  5.0 
IF  ( SUM9  -  180.0)  300,  3 
F0RMATI41H0N0RMALIZED  EX 
WRITE  OUTPUT  TAPE  4,316 
FORMAT  { 8 1  HO 
IMA  WBAR  WBARB 

WRITE  OUTPUT  TAPE  4,  131 

AO  I V  =  DIV*WFUNCTF (COEFT 
00  318  J  =  1,1 
AB1  =  B  1  l J ) / AD  I  V 
AB2  =  B2 1 J ) / AD  I  V 
ASBB  =  SBB ( J  )  /AD  I V 
FORMAT  I F20 . 2 ,  3F20.5/) 
WRITE  OUTPUT  TAPE4 ,  317, 
GO  TO  336 

WRITE  OUTPUT  TAPE  4,  296 
SUM  1 0  =  180.0 


FT  I  3 ) ,  090.0) 

ALPHA2 ,  ALPHA3 ,  SUM9  ) /WFUNCTF { ALPHA  1 

,  SUM9 ,  WFUNCT ,  VFUNCT 

00,  316 

PERI  MENTAL  WBARS  AND  SIGMAS/) 

ANGLE  WBAR 

/) 

6 

(1),  COEFT ( 2  )  ,  COEFT ( 3 ) ,  090.0  ) 


B3 { J ) ,  AB 1 ,  AB2 ,  ASBB 


t 


SIC 


324  FORMA  T  (  F  7 . 2  »  F10.6,  F20.6/) 

325  WFUNCT  =  WFUNCTF  (COEFTM  )  ,  COEFT(2),  COEFT13),  SUM  I  0 ) /WFUNCTF 
1 ( COEFT ( 1 ) *  COEF I ( 2 ) ,  COEFT(3),  270.0) 

335  VFUNCT  =  WFUNCTF { ALPHA  1 ,  ALPHA2,  ALPHA3,  SUM  1 0 ) /WFUNCTF ( ALPHA  1  , 

1  ALPHA2 »  ALPHA3 ,  270.0) 

WRITE  OUTPUT  TAPE  4,  324,  SUMIO,  WFUNCT,  VFUNCT 

SUM  1 0  a  SUMIO  +  5.0 

IF  (SUMIO  -  270.0)  325,  325,  635 

635  FORMAT! 4 IHONORMALI ZED  EXPERIMENTAL  WBARS  ANO  SIGMAS/) 

WRITE  OUTPUT  TAPE  4,635 

1635  FORMAT  (  8 1  HO  ANGLE  WQAR  S 

IMA  WBAR  WBARB  /) 

WRITE  OUTPUT  TAPE  4,  1635 

ADIV  =  DIV*WFUNCTF{COEFT{ I ) ,  COEFT ( 2 ) ,  COEFT { 3 ) ,  270.0) 

DO  637  J  =  1,L 
AB1  =  BKJJ/ADIV 
AB2  =  B2( J  J/ADIV 
ASBB  =  SBB ( J ) /ADI  V 

636  FORMAT ( F20 .2 ,  3F20.5/) 

637  WRITE  OUTPUT  TAPE4,  636,  B3(J),  AB 1 ,  AB2 ,  ASBB 

336  FORMAT { 2H0  //////) 

338  JK  =  2 

GO  TO  201  ; 

340  FORMAT { 52H0  ANISOTROPY-P  ANISOTROPY-C  SIGMA/) 


I  9H0SIGMA  =  ,  F6.4/) 

R V 1  =  ( ( COEFT( 2 )  +  COEFT (3) )/COEFT( l ) )**2 

PAREN1  =  (  SC0EFT(2)**2  +  SCOEFT(3)*»2  +  2. 0»CEE  I  NV  (  2 , 3)5/  DI V*  *2  )  / 
1  ( ( COEFT ( 2 )  +  C0EFT(3)) **2) 

PAREN2  =  ( SCOEFTI 1 )/COEFT( 1 ) )**2 

PAREN3  =  (-2.0/DIV**2)M  (CEEINVt  1 ,2)  +  CEE  I NV  (  2 , 3  )  )  /  (  OLEFT  (  1  )  * 

I  ( COEFT ( 2  )  t  COEFT ( 3 ) ))  ) 

SANISO  =  SQRTF(RV1*EPSSQ*(PAREN1  +  PAREN2  +  PAREN3 ) ) 

WRITE  OUTPUT  TAPE  4,  340 

WRITE  OUTPUT  TAPE  4,  241,  AN  I  SO (  1 ) »  ANIS012),  SANISO 

341  FORMAT (42H0C0EFFICIENTS  AND  SIGMAS  FOR  COSINE  SERIES/) 

WRITE  OUTPUT  TAPE  4,  341 

342  FORMAT(2F10.5/) 

Ui 

vO 

WRITE  OUTPUT  TAPE  4,  342,  (COEFT(I),  SCOEFT(l),  I  =  1,5) 

WRITE  OUTPUT  TAPE  4,  336 

RETURN 

END 

END 


FORTRAN  PROGRAM 


Appendix  Part  B 


Solid  Angle  Correction 


60 


■1 


PROGRAM  SOLANG 

409  FORMAT ( I5H0PR0GRAM  SOLANG/) 

410  FORMAT ( 82 HOTHE OR ETICAL  CORRELATION  FUNCTION  COEFFICIENT  CORRECTION 
IS  FOR  FINITE  SOLIO  ANGLE  /) 

411  FORMAT ( 2 1 HOMETHOD  OF  M.  E '.  ROSE/) 

WRITE  OUTPUT  TAPE  4,  409 

WRITE  OUTPUT  TAPE  4,  410 
WRITE  OUTPUT  TAPE  4,  411 

01  MENS  ION  H(2),T(2),R(2),BP(2) , GAMMA ( 2 ) , TAU ( 2 ) , A  INTO ( 2 ) , 

IAINT2I2)  ,AINT4{2) ,  DELTA ( 2 )  >  OELTAP (2 ) 

POLYOF(D)  =  D/D 

P0LY2F  {  0 )  =  C.5M3.0»(C0SF(D)**2)  -  1.0) 

P0LY4F ( D  )  =  0. 125»(35.0* (C0SF(D)**4)  -30 .0  * ( COSF ( D ) * *2 )  +3.0) 

EPSF { A  *  E  »  C  )  =  1.0  -  EXPF {-A*E*(  1 .O/COSFIC ) ) ) 

EPSPFt A,B,C,D)  =  1.0  -  E  XPF l —A* (B*(  l.O/SINFIC))— D*(  1. O/COSFIC) )) ) 
ARGOF { A  t  E  t  C )  =  POL YOF { C) *S INF ( C ) *EPSF ( A , E , C ) 

ARGPOF { A»  B  »C  »D  )  =  POL YOF ( C ) *SI NF ( C ) *EPSPF ( A , B , C , D ) 

ARG2F { A  » E  » C )  =  POL Y2F I C) * SINF ( C ) *EPSF { A , E , C ) 

ARGP2F ( A , B  » C  » D )  =  POL Y2F ( C ) * S I NF ( C ) *EPSPF ( A , B , C , D ) 

ARG4F ( A , E  »  C )  =  POL Y4F { C ) * S I NF ( C ) *E PSF ( A , E , C ) 

ARGP4F(A,B,C,D)  =  POL Y4F { C ) *S I NF I C ) *EPS PF I  A , B , C , 0 ) 

N  =  1 


210  FORMAT  (4F10.6) 

215  READ  INPUT  TAPE  3,  210,  R(1),  T(l),  H(l),  TAU(l) 

READ  INPUT  TAPE  3,  210,  R12),  T{2),  H(2),  TAU ( 2 ) 

220  DIV  =  200.0 
230  00  405  1=1,2 

BP ( I )  =  ATANF { R { I ) / ( H ( I )  +  Til))) 

GAMMA ( I )  =  ATANF { R { I ) /H( I ) ) 

DELTA (  I  )  =  BP { I  )  /D  I V 
DELTAP ( I )  =  ( GAMMA { I )  -  BP(I))/DIV 
240  BETA  =  0.0 

SUMO  =  ARGOFITAUII),  T(I),  0.0)  +  4 . 0* ARGO Ft  TAU ( I ) ,  T ( I )  ,  DELTA(I) 
1  ) 

250  SUMO  =  SUMO  +  2.0*ARGOF{ TAU( I ) ,  T(I),  (BETA  +2.0*DEITA{ I ) ) )  +  4 .0  * 

1 ARGOF ( TAU (  I ) ,  T ( I ) , ( BETA +  3. 0*DEITA ( I  )  )  ) 


BETA  =  BETA  +  2.0*DELTA(I) 

IF  ((BETA  +  3.0*DELTA( I) )  -  BP ( I ) ) 250 , 250 , 260 
260  SUMO  =  SUMO  +  ARGOF ( TAU( I ) ,  T( I ) ,  BP(I)) 

SUMPO  =  ARGPOF { TAU ( I ) ,  R(I),  BP(I),  H(I))  +4 . 0* ARGPO F ( T AU ( I ) , 
1 R (  I ) ,  BP ( I  )  +  DELTAP (  I  )  ,  H(I))  ; 

BETA  =  BP (  I  )  +  2.0*DELTAP( I ) 


270  SUMPO  =  SUMPO  +2.0*ARGP0F ( TAU( I ) ,  R(I),  BETA,  H(I))  +4.0*ARGP0F 
1 ( TAU { I ) ,  R ( I  )  ,  BETA  +DELTAP( I ) ,  H C I ) ) 

BETA  =  BETA  +  2.0*DELTAP(I) 

IF  {(BETA  +  DELTAP ( I ) )  -  GAMMA ( I )) 270 , 270 , 280 
280  SUMPO  =  SUMPO  +  ARGPOF (TAU( I ) ,  R(I),  GAMMA  II),  H  C I ) ) 

290  AINTO(I)  =  (DELTA (I)/3.0)*SUM0  +  ( DELTAP ( I ) /3 . 0 ) *SUMP0 
300  BETA  =  0.0 

SliM2  =  ARG2F  (  TAU  (  I  )  ,  T(I),  0.0)  +  4 . 0*  ARG2  F  (  T  AU  (  I  )  ,  T(I),  DELTA(I) 
1  ) 

3)0  SUM2  =  SUM2  +  2. 0*ARG2F{ TAU ( I ) ,  T(I),  (BETA  +2 . 0*DEL TA (I) ) ) +4 . 0* 
1ARG2F{ TAU( I ) ,  T( I ) , ( BETA +3 . 0*DELTA ( I ) ) ) 

BETA  =  BETA  +  2.0*DELTA(I) 

IF  ({BETA  +  3. 0*DEL  TA { I) )  -  BP ( I ) ) 3 1 0 , 3 1 0 , 320 
320  SUM2  =  SUM2  +  ARG2F ( TAU( I ) ,  T(I),  BP(I)) 

SUMP2  =  ARGP2F ( TAU ( I ) ,  R(I),  BP ( I ) ,  H ( I ) )  +4 . 0« ARGP2 F ( TAU (I) , 

1 R ( I ) ,  BP ( I ) +DELTAP ( I ) ,  H(I)) 

BETA  =  BP ( I  )  +  2.0*DELTAP ( I ) 

330  SUMP2  =  SUMP2  +2.0*ARGP2F(TAU(I ) ,  R(I),  BETA,  H ( I ) )  +4.0*ARGP2F 
1 l TAU ( I ) ,  R (  I  )  ,  BETA  +DELTAP ( I ) ,  H(I)) 

BETA  =  BETA  +  2- 0*DELTAP ( I  ) 

IF  {(BETA  +  DELTAP { I ) )  -  GAMMA ( I  )) 330 , 330 , 340 
340  SUMP2  =  SUMP2  +  ARGP2F (TAU ( L ) ,  R(I),  GAMMA ( I ) ,  H ( I ) ) 

350  A I NT2 ( I )  =  ( DELTA( I ) /3.0 )*SUM2  +  ( DELTAP ( I ) /3. 0 ) *SUMP2 


360  BETA  =  0.0 

SUM4  =  ARG4F ( TAU ( I ) ,  TCI),  0.0)  +  4 . 0*ARG4F l TAU (  I ) ,  TCI),  CELT  A{  I ) 
1  ) 

370  SUM4  =  SUM4  +  2. 0*ARG4F( TAU ( I ) ,  TCI),  (BETA  +2 . 0 »DELTA (I) ) ) +4 . 0* 

1 ARG4F ( TAU ( I  )  ,  TCI) , ( BETA +3 . 0*DEL TA ( I ) ) ) 

BETA  =  BETA  +  2.0*DELTA(I) 

IF  ((BETA  +  3.0*DELTA( I) )  -  BP ( I )) 370, 370 , 380 
3£0  SUM4  =  SUM4  +  ARG4F  (  TAU(  I  )  ,  T  (  I  )  ,  BP(D) 


SUMP4  =  ARGP4F ( TAU ( I ) ,  R(I),  BP ( I ) ,  H ( I ) )  +4 . 0 *ARGP4 F ( T AU ( I ) , 

IRC  I ) ,  BP ( I ) +DELTAP ( I ) ,  H ( I ) ) 

BETA  =  BP (  I  )  +  2 . 0*DELTA  P ( I ) 

390  SUMP4  =  SUMP4  +2 . 0 * ARGP4 F ( TAU (I) ,  R(I),  BETA,  H( I ) )  +4.0»ARGP4F 
UTAU(I),  R  (  I  )  ,  BETA  +DELTAPU),  H(I)) 

BETA  =  BETA  +  2 . 0#DELTAP ( I ) 

.  IF  ((BETA  +  DEL  TAP ( I ) )  -  GAMMA { I )) 390 , 390 , 400 
400  SUMP4  =  SUMP4  +  ARGP4F ( T AU ( I ) ,  R ( I ) ,  GAMMA ( I ) ,  H ( I ) ) 

405  A  I NT4 (  I  )  =  (DELTA! I)/3.0)*SUM4  +  (DELTAPt I )/3.0)*SUMP4 
AN0RM=  1.0/(AINT0(  1 1  *A INTO  C  2  >  > 

AQ2  =ANORM*A INT2 (  l)*AINT2(2) 

AQ4  =AN0RM*AINT4(  1  )*AINT4(2) 

412  FORMAT!  17H0INPUT  PARAMETERS/) 


413  FORMAT! 89H0CRYSTAL  NUMBER 


RAO  I  US 


THICKNESS 


DISTANCE  TO  SO 


1URCE  LINEAR  ABSORPTION  COEFF  /) 

414  FORMAT (19,  F16.6,  F13.6,  F17.6,  F21.6/) 

4 1 5  FORMAT (8H0RE  SUL  TS/ ) 

420  FORMAT ( 6H0Q2  =  F9.6/) 

430  FORMAT ( 6H0Q4  =  F9.6/) 

440  FORMAT! 30H0ENO  OF  PROGRAM  SOLANG  RESULTS/) 

WRITE  OUTPUT  TAPE  4,  412 

WRITE  OUTPUT  TAPE  4,  413 

WRITE  OUTPUT  TAPE  4,  414,  (I,  R(I),  TII),  HID,  TAUII),  l  =  1,2) 

WRITE  OUTPUT  TAPE  4,  415 

WRITE  OUTPUT  TAPE  4,  420,  AQ2 

445  WRITE  OUTPUT  TAPE  4,  430,  AQ4 

N  =  N+  1 

IF  (N-17)  215,215,520 
520  WRITE  OUTPUT  TAPE  4,  440 
END  FILE  4 
PAUSE 
END 
END 
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PROGRAM  CRT  I  ME  „  \ 

DIMENSIONAN U 500 ) , AN2( 500) ,  CT(500),  TAUDK500),  TAUD2(500),  2 

1  T ( 500 ) t  T AUC ( 500  )  ,  STAUC ( 500 ) , NRUM ( 500 )  3 

398  FORMAT (23H026  APR  CRT  I  ME  MOO  > 

WRITE  OUTPUT  TAPE  4,  398  ^ 

40  FORMAT  (13)  3-1 

45  READ  40,  MM  3-2 

LM  =  l  3-3 

50  FORMAT  (13)  ** 

READ  50,  L 

60  FORMAT ( 3F 12. 1 ,  2F12.9,  FIO.l,  110)  3 

READ  60, (AN1 (K) ,AN2(K) ,  CT(K),  TAUDI(K),  TAUD2 ( K ) ,  T ( K ) , NRUN ( K ) , K=6 

11, L) 


DO  70  K=  I  , L 
AND0T1  =  AN  1  ( K ) /T ( K ) 

AND0T1P  =  ANDOT  1  * (  1.0  +  TAUO  l  ( K ) *ANUOT 1 ) 
AND0T2  =  AN2 ( K ) / T ( K ) 

AND0T2P  =  AND0T2* (1.0  +  TAUu2( K ) *AND0T2 ) 
TAUC(K)  =  (CT(K)/T(K))/(2.0*ANDOTIP*ANDOT2P) 

70  STAUC(K)  =  TAUC(K)/SQRTF(CT(K) ) 

SUM  1  =  0.0 

DO  71  K  =  1,  L 

71  SUMI  =  SUM  1  +  TAUC(K)/(STAUC(K)**2) 

SUM2  =*0.0 

DO  72  K  =  1,  L 


8 

9 

1U 
1  1 
12 

13 

14 

5 

16 


o 

oo 


72  SUM2  =  SUM2  +  1.0/(STAUC(K)**2) 

T  AUCM=  SUM  1 / SUM2 
D  =  L 

STAUCM  =  SQRTF (  1  . 0/SUM2 ) 

90  FORMAT ( 29H0  COINCIDENCE  RESOLVING  TIMES/) 

WRITE  OUTPUT  TAPE  4,  90 

100  FORMAT ( 3 1H0TAU  AND  SIGMA  TAU  FOR  EACH  RUN/) 

WRITE  OUTPUT  TAPE  4,  100 

110  .FORMAT  (9P2F20.5,  110) 

WRITE  OUTPUT  TAPE  4,  I  1 0  » ( TAUC ( K ) ,  S TAUC ( K ) , NRUN { K ) , 
120  FORMAT {57H0MEAN  RESOLVING  TIME  SID  DEV 

_ 1 R / 1  _  _  _  ... 

SUM3  =  0.0 
DO  130  K  =  1,  L 

130  SUM3  =  SUM 3  +  ( (TAUC(K)-TAUCM)/STAUC (K)  )»*2 
RMS  =  SQRTF ( SUM3/SUM2) 

WRITE  OUTPUT  TAPE  4,  120 

140  FORMAT  (9P3F20.5) 

WRITE  OUTPUT  TAP  E4  »  140,  TAUCM , STAUC M , RMS 

150  LM  =  LM  +  1 

IFILM  -  MM)  50,  50,  160 

160  END  FILE  4 
PAUSE 
END 
END 


1 

1* 

2S 

21 

22 

23 
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23 
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27.  1 
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33 
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..JOB  PHELPS 

GAMMA-GAMMA  ANGULAR  CORRELATION  ANALYSIS 
10  MAY  VAA  11/18/ 

NUMBER  OF  SETS  OF  DATA  TO  BE  ANALYZED 

2 

DATA  SET  NUMBER 


CHECK  INPUT  PARAMETERS 

TAUC  STAUC  TAUD1  TAUD2  BACK  1  BACK2  T 

.00000025786  .0000000028  .0000035000  .0000035000  3.8000  3.5700  7200.00 

o 

THEORETICAL  LEGENDRE  POLYNOMIAL  COEFFICIENTS  CORRECTED  FOR  SOLID  ANGLE 
l. 000000  .099615  .008368 

NUMBER  OF  RUNS  AT  EACH  ANGLE 
90.00  3 

105.00  3 

120.00  .  3 
135.00  3 

150.00  3 

165.00  3  (, 

180.00  3 


l 


CHECK  DATA  INPUT  OF  Nl,  N2,  AND  CT  AND  INPUT  PARAMETERS 


ANGLE- 


90.00 


N  J 

7715372.0 

7695190.0 

7676042.0 


N2 


CT 


TAUC 


STAUC 


TAUD1 


TAU02 


BACK  1 


BACK2 


8493912.0 

8483583.0 

8421415.0 


10525.0  .00000025786  .0000000028  .0000035000  .0000035000  3.300C  3.5700  7200.00 
10214.0  .00000025786  .0000000028  .0000035000  .0000035000  3.3000.  3.5700  7200.00 
10118.0  .00000025786  .0000000028  .0000035000  .0000035000  3.3000  3.5700  7200.00 


ANGLE-  105.00 


Nl 

7712620.0 

7735075.0 

7759638.0 


N2 


CT 


TAUC 


STAUC 


TAUD1 


TAUD2 


BAIX  1 


8520243.0 

8543918.0 

8556720.0 


10344.0  .00000025786  .0000000028  .0000035000  .0000035000  3.3000 
10432.0  .00000025786  .0000000028  .0000035000  .0O0008S00U  3.3000 
10485.0  .00000025786  .0000000028  .0000035000  .0000035000  3.3000 


BACK2 

3.5700 

3.5700 

3.5700 


ANGLE-  120.00 


Nl 

7696778.0 

7701444.0 

7744192.0 


N2 


8545595.0 

8523666.0 

8534320.0 


CT 

10420.0 
10422.0 
1057  1.0 


TAUC 

.00000025786 

.00000025786 

.00000025786 


STAUC 

.0000000028 

.0000000028 

.0000000028 


TAUD  1 

.0000035000 

.0000035000 

.0000035000 


TAUD2 

.0000035000 

.0000035000 

.0000035000 


BACK  1  BACK2 


3.3000 

3.3000 

3.3000 


3.5700 

3.5700 

3.5700 


7200.00 
7200. OC 
7200.00 


7200.00 

7200.00 

7200.00 


ANGLE-  135.00 


N  1 

771  1775.0 
7796272.0 
7760244.0 


N2 

8522919.0 

8552850.0 

8476438.0 


CT  TAUC 

10745.0  .00000025786 
10781.0  .00000025786 
10822.0  .00000025786 


STAUC 

.0000000028 

.0000000028 

.0000000028 


T  AUD  1 

.0000035000 

.0000035000 

.0000035000 


TAUD2  BACKFi 

.0000035000  3.8000 
.0000035000  3.8Q30 
.0000035000  3.8000 


BACK2  T 

3.5700  7200.00 
3.5700  7200.00 
3^5700  7200.00 


ANGLE-  150.00 


N  1 

7665166.0 

7723064.0 

7551069.0 


N2 

8540160.0 

8581780.0 

8418299.0 


CT  TAUC 

10734.0  .00000025786 
11152.0  .00000025786 
10557.0  .00000025786 


STAUC 

.0000000028 

.0000000028 

.0000000028 


T  AUD  1 

.0000035000 

.0000035000 

.0000035000 


TAUD2  BACK! 

.0000035000  3.8 DJ 0 
.0000035000  3.8C0G 
.0000035000  3.8030 


BACK2  T 

3.5700  7200.00 
3.5700  7200.00 
3.5700  7200.00 


ANGLE-  165.00 


N  1 

7655554.0 

7726060.0 

7663836.0 


N2 

8519721 .0 
8546518.0 
8565178.0 


CT  TAUC 

11105.0  .00000025786 
11214.0  .00000025786 
11203.0  .00000025786 


STAUC 

.0000000028 

.0000000028 

.0000000028 


T  AUD  1 

.0000035000 

.0000035000 

.0000035000 


TAUD2  BACK! 

.0000035000  3.8050 
.0000035000  3.8010 
.0000035000  3.8030 


BACK2  T  . 

3.5700  7200.00 
3.5700  7200.00 
3.5700  7200. UU 


ANGLE-  180.00 


N  1 

7716070.0 

7724981.0 

7713354.0 


N2 

8530187.0 

8524972.0 

8583727.0 


CT  TAUC 

11177.0  .00000025786 
11083.0  .00000025786 
11206.0  .00000025786 


STAUC 

.0000000028 

.0000000028 

.0000000028 


T  AUD  1 

.0000035000 

.0000035000 

.0000035000 


TAU02  BACKS 

.0000035000  3.8010 
.0  0  0  0  0  3  5  0  0  0  3.8  030 
.0000035000  3.8030 


BACK2 

3.5700 

3.5700 

3.5700 


T 

7200.00 

7200.00 

7200.00 


1 


RESULTS  OF  ANALYSIS 


NDOT 1  PRIME 

ND0T2PRIME 

CTOTALOOT 

ANGLE-  90.00 

1075.5984 

1  184.5810 

1  .46  18 

1072.7744 

1 183. 1346 

1 .4186 

1070.0951 

1174.4292 

1 .4053 

ANGLE-  105.00 

1075.2133 

l  188.2683 

1.4367 

1078.3555 

1  191.5838 

1.4489 

1081. 7927 

1 193.3766 

1 .4562 

ANGLE-  120.00 

1072.9966 

1191.8187 

1.4472 

,1073.6495 

1 188.7477 

1.4475 

>1079.6313 

1 190.2397 

1.4682 

ANGLE-  135.00 

1075.0951 

1 188.6431 

1 .4924 

1086.9193 

1  192.8347 

1.4974 

1081.8775 

1 182. 1340 

-  1.5031 

ANGLE-  150.00 

1068.5732 

1  191.0575 

1.4908 

1076.6748 

1 196.8862 

1.5489 

1052.6092 

1 173.9929 

1  .'4662 

CACCDOT 

CGENOOT 

W 

SIGN 

.6571 

.8047 

.63573E— 06  - 

.14678 

.6546 

.7640 

.  60594E— G’6 

.14590 

.6481 

.7571 

.60645E— 06 

.14654 

.6589 

.7778 

.  6  1 275E— 06 

.14575 

.6627 

.7862 

.6  1587E—C6 

.14551 

.  6658 

.7905 

.6  16306-06 

.14523 

.6595 

.7877 

.6200  16- C  6 

.  14596 

.6582 

.7893 

.622496-06 

.14618 

.6627 

.8055  . 

•  63093E— 06 

.14607 

.6590 

.8333 

.  65638E— 06 

. 14736 

.6686 

.8287 

.6433  6E-  C  6 

. 14597 

.6596 

.8435 

.663866-06 

.14759 

.6564 

.8345 

.65996E-06 

.14775 

.6646 

.8843 

.6907  IE-0  6 

.  14810 

.6373 

.8289 

.675286-06 

.15022 

ANGLE-  165.00 


1067.2283 

1077.0940 

1068.3871 


1188. 1952 
l 191.9479 
1 194.5611 


1.5424 

1.5575 

1.5560 


.6540 
.662  l 
.6582 


.8884 

.8954 

.8978 


.  70520$.— 06 
.  7020  IE-06 
.  7  0808E-C6 


.14966 
.1487 
.  1493 


ANGLE-  180.00 

1075.6961 

1076.9430 

1075.3161 

1 189.6609 
1188.9306 
1  197.  1589 

1.5524 

1.5393 

1.5564 

.6600 

.6603 

.6639 

T  8924 
.8790 
.8925 

70  1 9  IE— 06 
69098E-06 
69783E-06 

ANGLE 

WBAR 

SIGMA  WBAR 

WDARB 

90.00 

•  6 l 599E-06 

. 84528E-08 

.98458E— 08 

105.00 

. 6 1 498E-06 

. 84003E-08 

.11  189E-08 

120.00 

. 62448E-06 

.84334E-08 

• 33054E-08 

135.00 

• 65445E-06 

. 84853E-08 

.59896E-08 

150.00 

. 67530E-06 

. 85840E-08 

•88763E-08 

165.00 

• 70509E-06 

.86  172E-08 

.  17529E-08 

180.00 

. 69690E-06 

.8581 2E-08 

.3 l 886E-08 

LEAST  SQUARE  FIT  COEFFICIENTS  OF  LEGENDRE  POLYNOMIALS  AND  THEIR  SIGMAS 
1.00000  .00565 

.09200  .01024 

.00991  .01239 

EPSILON  SQUARE  =  .644725E+00 


.1489 
.  1485 
.1484 


roNitri  inift 


EXPERIMENTAL  FUNCTION 


90.00 

1.000000 

95.00 

1.000802 

100.00 

1.003216 

105.00 

1.007256 

110.00 

1.012936 

115.00 

1.020250 

120.00 

1.029152 

125.00 

1.039539 

130.00 

1.051232 

135.00 

1.063963 

140.00 

1.077377 

145.00 

1.091035 

150.00 

1.104433 

155.00 

1.1 17029 

160.00 

1.128273 

165.00. 

1. 137647 

170.00 

1.144698 

175.00 

1.149077 

180.00 

1.150562 

THEORETICAL  FUNCTION 
I .000000 
1.000943 
1.003769 
1.008467 
1 .015010 
1.023340 
1.033355 
1.044892 
1.057716 
1 .07  151 1 
1 .08S886 
1 .  100376 
1 . 114467 
1.127615 
1 . 139280 
1 . 148956 
1 .156209 
1 . 160702 
1  .  162223 


NORMALIZED  EXPERIMENTAL  WBARS  AND  SIGMAS 


ANGLE 

90.00 

105.00 

120.00 

135.00 

150.00 

165.00 

180.00 


WBAR 
1.00623 
1.00458 
1.02010 
1.06907 
1.  1031  1 
1.15178 
1. 13840 


ANISOTROPY-P  ANISOTROPY-C 

.150562200670  .150562197679 


COEFFICIENTS  AND  SIGMAS  FOR 

l. 00000  .01053 

.10529  .05816 

.04527  .05658 


COSINE  SERIES 


I GMA  WBAR 
.01381 
.01372 
.01378 
.01386 
.01402 
.01408 
.01402 


W8ARB 
.  0 1608 
.00183 
.00540 
.00978 
.01450 

.002  m 

.00  52T 


SIGMA 

.028830656297 


Appendix  E  FORTRAN  Program 
Glossary  of  Terms  for  the 
Least  Squares  Analysis 

1.  Check  of  input  parameters 

TAUC  Coincidence  unit  resolving  time  (sec) 

STAUC  Standard  deviation  of  T  (sec) 

TAUDl  System  dead  time,  Channel  1  (sec) 

TAUD2  System  dead  time,  Channel  2  (sec) 

BACK1  Background,  channel  1  (counts  per  sec) 

BACK2  Background,  channel  2  (counts  per  sec) 

T  Duration  of  run  (sec) 

20  Check  data  input  of  N^N^and  CT  and  input  parameters 

Total  counts,  channel  1 

Total  counts,  channel  2 

CT  Total  coincidence  counts 


3.  Results  of  analysis 

ND0T1PRIME  Count  rate,  channel  1,  corrected  for  system  dead 
time  and  background  (counts  per  sec) 

Corresponds  to  ft.,  in  text 
ND0T2PRIME  Same,  channel  2  1 


CTOTALDOT  CT 

CACCDOT  C A 
A 


CGENDOT 

W 


SIGMAW 


W  (0.)  - 

ex'  i 


*1*2 


/  i 


W  (G.)l  = 

ex  i' 

X 

w  (@.) 

ex  iy 

c 

L  J 

G 

- 

WBAR 


W_  (&:)- 

e  x 


I 

1  IkS* 


a 


where  N  is  the  number  of  runs  at  the  angle  0^. 

EXPERIMENTAL  FUNCTION  W.  (9)  normalized  to  Wn  (90)  and 

Is  Is 

evaluated  at  5°  intervals. 

THEORETICAL  FUNCTION  Wth(Q)  normalized  to  W  (90)  and 

evaluated  at  5°  intervals. 


NORMALIZED  EXPERIMENTAL  WBARS  AND  SIGMAS 

WBAR  W  (0.)  normalized  to  W  (90) 
ex  i  J.  s 


ANISOTROPY  -  P 


W.  (180) 


1 


based  upon  Legendre  series 
representation  of 


ANISOTROPY  -  C 


Same 
of  W 


,  based  on  Cosine  series  representation 
Is  ’ 


SIGMA  Standard  deviation  of  ANISOTROPY  -  P 
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